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ABSTRACT 


Performance  measures  and  associated  criteria,  for  linear  constant  coef¬ 
ficient  system  forced  by  deterministic  input?  are  investigated,  with  particu¬ 
lar  reference  to  flight  control  system.  I*  in  shown  that  the  application  of 
performance  treasures  is  facilitated  by  substituting  for  the  actual  flight  con¬ 
trol  system  an  equivalent  lev- order  linearised  systes  having  stellar  dynamic 
characteristics.  A  critical  survey  of  current  performance  neasurcs  is  given, 
and  nev  methods  for  the  analytic  evaluation  of  socc  irxlic ini  error  measures 
are  prerented. 

Hunercus  criteria  are  examined  with  regard  to  their  validity,  selectivity, 
and  case  of  application,  Korralized  presentations  are  used  so  that  practical 
limitations  on  the  tine  scale  of  the  response  (c.g.,  due  to  pover/inertia  re¬ 
strictions)  can  be  taken  into  account  separately.  It  is  concluded  that  nininua 
ITAE  (integrated  ticc  moment  of  absolute  error)  and  minima:  ITE2  (integrated 
time  accent  of  error-squared)  have  particular  merit.  The  ITAE  criterion  yields 
smooth  indicial  responses  having  little  oversnoot,  but  its  analytic  description 
is  complicated.  Of  the  criteria— minims*  IE2  (integrated  error-Jfjiarcd) ,  mini¬ 
mum  ITE2,  IT^E2  (integrated  first-,  second-,  and  third  time  accents  of 

error-squared)  — IE2  bar.  simple  analytic  forrs,  but  selects  poor  indici-l  re¬ 
sponses;  IT5e2  responses  are  as  good  as  those  selected  by  ITAE,  but  IT^E2  (and 
also  USE2)  aj»al/tlc  expressions  arc  too  cccplicatcd  for  general  use.  ITE^  se¬ 
lects  moderately  smooth  and  veil -damped  responses  (less  good  than  ITAE),  but 
possesses  tractable  analytic  fores.  Therefore,  ITE2  is  recccnendcd  for  analytic 
investigations,  whereas  ITAE  is  prtfen-d  for  optimization  using  analog  computers. 
Sooc  other  criteria  also  appear  promising  for  us*'  in  conjunction  with  digital 
cocputers,  but  require  further  investigation  to  determine  t'ueii  validity  and 
selectivity. 
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A(t)  «  (ip  Appendix  B) 

*0  *»  2ke"°^l/K  0n  Appendix  B) 

[a]  MatrJx  of  coefficients  (cee  Tab)*  T) 
a  Coefficient  of  aperiodic  ters  In  er^>r  reaponae 

•0  Arguxent  of  Gaaea  function  (»  »/?) 

aj  Argument  of  G*eo%  luwtion  (*  5/2) 


a?  Argusx-nt  of  Oaara  function  (•  ’>/2) 
aj  Coefficient  of  ■*  in  characteristic  equation,  &(•)  •  o 
ap  Coefficient  of  an“*  ip  denominator  cf  fc(t) 


3  Set  Eq  E-57  and  B-  JO  (Appendix  B) 

B{r,ao>  »  l'(s)r(g(j)/r(*  ♦  pq),  the  Beta  function 
b  -  ee"gt 

b  Coefficient  of  c*  in  Eq  A? 

bQ  Koran 11 zed  coefficient  of  c*  in  denominator  of  C(»)/R(e) 


C(a) 

c(t.o 

C(l,a) 

?! 


-  c*8T 

- 


See  Eq  BO  8  (Appendix  B) 

-  c(t,;)sc 

»a-Joo 


P»th»  of  tatogrttloo  (u*  rig.  B-}) 


c  .  y/TT£  (in  Appendix  B) 

c  Coefficient  of  a  la  Eq  A7 


c0  Bon*%li»ed  c<rriot...t  rf  -  in  dcn-slnator  of  C(n)/R(a) 
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c;  ■  (la  Appendix  E) 

Cp  Peak  overshoot 

c(t)  Kesperse 

cc(t)  Response  envelope 

D  ?**ancfer  function  denccticator 

Coefficient  cf  cn“*  in  niraerator  of  K(s) 

s(*>  *£Wtfl 

Za(») 

tj.,  1th  error  coefficient 

e  -  2.718 

s(t>  Error  response 

|e(t)|  Absolute  \alue  of  e(t) 

u(«) 

f-(®)  -,£[*-(*)] 

C(s)  Transfer  function 

Cxr  -  C(9)/R(*)  -  C(t)/l  C(«) 

Cre  -  E(c)/R(s)  -  1/1  *  GU) 

jG_|  Deviation  ratio 


5 

50 

Ho 


b(x) 


h. 


IAE 


IE 


*  7  -  a 

ArgWwrot  of  Beta  function 
*Jcc  Eg  B-U  (Appendix  B) 

See  Eg  B-83  (Appendix  3) 

■  h+  «  h.,  the  transfer  characterise* c  cf  the  full-vave  linear  .-ectificr 

That  part  of  transfer  characteristic  for  x  positive  or  negative, 
respectively 

Integral  of  absolute  error.  ^jf®  (e't)]dt 
Integral  of  trror,  e(t)dt 
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IE2 

ITAE 


Integra),  of  error- squared,  jT®  e(t)  2Ct 

Integral  of  time -weighted  absolute  error,  JP®  tje(c)  |dt 

Integml  of  tlae- weighted  error,  te(t)dt 

Integral  of  time-weighted  error- squared,  t  e(t)  *dt 

Integral  of  tt'  e- squared  times  error-squared,  t*  r(t)  *4t 

Integral  cf  tlae-cubed  .Lars  error- squared.  r.’  c.tj  'ct 

Integral  of  0th  tlneooatnt  of  error,  ^fl0 

•C°  Mt 

t°U!t 

Integral  of  a  function  of  error  not  involving  Ume  explicitly,  Uut 
Modified  Bessel  function  of  the  first  kind 
Bescel  function 


•V* 


See  Eq  B-19  (Appenlix  *) 


Coefficient  of  oscillatory  tsr*  in  error  response  of  trum-order  system 
(Appendix  B) 

Sea  Eq  B-61  (Appendix  B) 

Laplace  transform  operator 

Inverse  Xaplace  transforn  operator 

honing  mcnent  derivative  with  respect  to  sideslip 

See  Eq  B-6l  (Appendix  B) 

Peak  ratification  ratio 

Cider  of  "taraorlct"  in  output  signal 
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If  Transfer  function  r.userotor 

R  Susber  of  cloc)r-*icc  cnci»'elcx*iits  of  -5  by  C(s) 

%  Yavlcg  Docent  ae.Jrative  with  respect  to  sideslip 

K„  Sec  Eq  B-Sj  (Appendix  B) 

n  Sussation  or  infinite  product  Index  for  Beta  function 

n  Sy.-tea  order 

f  Iteber  of  poles  or  0(s)  in  rl*ht- half  plane 

?i  Coefficient  of  X*  lb  the  numerator  of  Eq  19 

p  Laplace  transfers  variable  if  1  is  independent  variable 

pj  Coefficient  of  X1  in  the  n  iterator  of  Eq  2} 

Qi  Ooefficien*  of  X*  in  the  dencalaator  of  Eq  19 

q  ExpoUfnt  (Appendix  B) 

qi  coefficient  of  X1  in  the  dennaiaato r  of  Eq  2J 

qA  Coefficient  of  s*  in  Cia)/R(s)  denominator 

*(•>  ^><0 

r  Suasatiot.  index 

r(t)  input 


s  Laplace  transfora  variable 

R  ■  sin”1  b 

*\  -  sin”1  bj 

T  A  tieje  constant 

Tc-  A  tis-  constant 

Tj  •[lc{aA1^7  -  o.  tiac  when  ratio  or  envelopes  equals  one 

Elevator  se. vc  tine  constant 
Teg  Equivalent  ti sx  coo»u*ni 
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r 


t 

*p 

IV 

u 

C 

V(tJ 

V 

w  or 
W(t) 

HaW 

>t(-) 


lo 

X 

x  or 
x(t) 

Y 

y 

y  or 
y(t) 


Closed-loo?  tine  constant 

Tlae,  or  ccK*?l2ed  tlxe 

•ficc  or  n^rrslired  tlae 

Delay  tlae 

Tlae  to  peak 

Rise  tice 

Settling  tlae 

See  Eq  B-fcfc  (Apj*:r.iix  Z) 

*  Re  (uj) 

An  error  function,  not  Involving  ti*e  explicitly 
-ke^ 

*  I*(®} 

See  Eq  B-M  (Appendix  B) 

jC 


S~-  In  B-69  (Apprallx  B) 


-  Be  U) 

Input  to  full-vave  linear  rectifier 
Open -loop  truisfet  function 
•  la  (s) 

Output  of  the  full-vave  linear  rectifier 
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Z  or 
=(t> 


Number  of  zeros  of  I  *  0(c)  in  rl(rht-balf  plane 
See  Sq  3-J7  (Appendix  S) 


a(s) 

3 


»(s) 

r(x) 

r 

> 

Ji  1 

^  I 

A. 


o 

5 

ie 

C- 

*0 

5 

; 

?i 

i 


**  x  ♦  ji  cctrplcx  variable 

Real  port  of  copies  roots  of  third -order  systen 
Jruscrator  of  C(s.i 

Ica^ir-ary  part  of  cccplex  rcotf  of  third-order  ayaten 
‘  ‘wF  -  t2  ■  of  (Appendix  A  only)  — i*'  u 
rvocainator  or  G(s) 

Cassa  function 

Real  root  of  third-order  3ystes  transfer  function  dencnlmtor 
laplaoe  tmnefors  variable  (in  Chapter  III) 

Characteristic  equation 


See  Eq  B-6'.  (Appendix  B) 

■  t*3'1  ~T^S 

c?  -  62 

■  -.an-1  2 

a 

•  <V(t) 

A  tisc  constant 
Elevator  deflection 

■  2  o  •  1 ,  ?,  ...  Hc'xtarn  factor 

»  l 


Seal  number  tlaa  -ere 

.  a XU)  •  i  •»  ji),  corpiex  variable 
Baxplng  ratio 
Closed-locp  dating  ratio 
•  ia  (?) 


nr 
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9  Angle  of  pitch 

8{t)  -  pt  ♦  ♦  -  | 

K  Open-loop  gain 

kc  Equalization  gain 

•(g  Gain  of  cervcoctcr  plus  amplifier  cucblnatlon 

k^/k  Gala  anrgln 

K^fK  Generalized  gain  aargia 

X  Open- loop  tiae  constant 

X  ■  %  ♦  o 

X  -  Oq* 

m  Magnitude  of  laplace  transfora  variable,  a 

Mc  Generalized  crossover  frequency 

t  *  V*  (O 

-l  Carping  ratio  of  closed-loop  roots 

o  Real  part  of  cc*>lex  variable,  s 

0|  A  particular  value  of  a 

t  Tiae  delay 

t  ■  t  -  T| 

0  •  ♦  ♦  6T,  (phase  angle) 

♦a  Phase  aargln 

0^  Generalized  phase  eargio 

-  tan*1  |  -  tan'1  y-^~a  (phase  angle) 
Xoimllzed  freqswtxy 

Taaginery  part  of  coeplcx  variable,  e  >  a  ♦  Joa 


t 

t 

°0 
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cd  Laplace  transfers  variable  vh*o  x  la  the  independent  variable 

a^j  mndvidih 

a^.  Croaaover  frequency 

aij  Closed- loop  root  natural  frequency 

a^(£)  See  £q  h-J2  (Appendix  B) 

cu  Undarjwd  natural  frequency 

<Dp  Peak  Magnification  frequency 

tq,  Frequency  of  instability 

4  Angle 

|  )  Absolute  value 

}  Magnitude  in  decibels 

d  Partial  derivative 

53  Susratloo 

II  product 
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Dynamic  performance  is  only  one  of  many  factors  to  1*6  considered  in 
assessing  the  merit  of  a  flight  control  ftjrsten-  Cost,  weight,  reliability, 
schedule,  etc.,  mat  alao  be  taken  into  account.  The  beet  choice  for  any 
given  requirement  can  be  found  only  by  weighting  each  of  these  factors  accord¬ 
ing  to  their  relative  importance.  Coat,  weight,  etc.,  are  measured  directly 
and  unequivocally  in  terns  of  dollars,  pounds,  etc.,  hut  at  the  present  time 
assessment  of  dynamic  performance  depends  heavily  on  intuition.  This  la  not 
due  to  any  shortage  of  performance  measures;  many  have  been  proposed.  The 
problem  la  which  measure  or  combination  of  measures  to  choose,  and,  having 
cade  the  choice,  hov  to  apply  the  resulting  criteria  to  the  system  under  con¬ 
sideration.  It  is  to  the  solution  of  this  problem  that  this  report  is  directed. 

The  work  reported  here  was  performed  under  an  Air  Force  contract  directed 
at  providing 

1 .  A  foundation  for  the  specification  of  dynviic  performance 
criteria  for  autcc&tle  flight  control  systems 

2.  The  methods  of  analysis  required  to  apply  such  criteria- 

It  was  convenient  tc  present  the  results  of  the  study  in  tve  parte.  The 
present  report  deals  with  perfammnee  criteria  appropriate  to  deterministic 
inputs;  rarri.cn  inputs  and  associated  topics  are  discussed  in  a  subsequent 
report. 

Chapter  I  consists  of  a  broad  discussion  of  dynamic  performance,  and  out¬ 
lines  exact  and  approximate  calculation  procedures.  It  Is  shown  thet  although 
actual  flight  control  systems  are  described  by  differential  equations  of  high 
order,  more  tractable  equivalent  systems  of  low  order  can  provide  a  convenient 
and  sufficiently  accurate  basis  for  analytical  optimization  techniques.  Quali¬ 
ties  defining  the  merit  of  the  dynamic  performance  of  given  systems  are  out¬ 
lined,  and  their  interpretation  into  numerical  performance  meaaures  is  dis¬ 
cussed.  Performance  criteria  are  then  defined  in  terms  of  optimal  values  of 
tnece  measures,  and  the  requirements  of  validity,  selectivity,  and  ease  of 
application  theft  a  gool  criterion  should  possess  are  formulated. 

Manuscript  released  by  authors  1  k  December  1961  for  publication  as  as  ASD 
Technical  Beport. 
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Chapters  II  and  171  survey  published  oeasures  or  dynsaic  performance. 
Cnapter  II  deals  primarily  with  criteria  rot  explicitly  invoicing  integral 
functions  of  error  (c-g.,  phase  oergin,  bandwidth,  rite  tine,  etc.),  and  pre¬ 
sents  extensive  correlations  of  tnese  criteria  for  second-order  systems. 

Chapter  III  is  concer^d  witn  indicinl  error  erasures;  i.e- ,  criteria  directly 
oeasuring  soeae  integrated  function  of  the  error  response  to  a  step  input.  This 
class  includes  sose  of  the  most  useful  criteria,  such  as  ITAE  (integrated  tine 
onr*.  of  absolute  error),  IE^  (integrated  error-squared),  etc.  Calculation 
procedures  for  these  criteria  are  presented;  in  particular,  analytical  expres¬ 
sions  for  ITAE  tvnd  soee  other  aeasures  not  previously  expressed  in  analytical 
*  Mto  -.r~'  published  for  the  first  tine.  Chapter  III  concludes  with  a  survey 
of  tlie  closel-loop  and  open -loop  characteristics  of  those  opt  Inal  systeta  which 
ninisisc  particular  indicia!  error  treasures. 

Copter  IV  derives  ge”«ral  exact  fora-lac  expressing  the  effect  of  a  pure 
tine  lag  on  system  i'idicial  error  measures.  By  combining  these  formulae  with 
tl*  observation  that  the  responses  of  high-order  good  ar.d  optical  systems 
closely  approximate  to  the  responses  of  lover-order  systems  with  a  tine  delay, 
a  metnod  of  approximating  to  performance  measures  of  high-order  systems  is 
devised.  Ar  exasple  is  given  to  show  how  this  approximation  simplifies  opti¬ 
misation  procedures. 

Chapter  V  generalizes  some  of  the  formulae  for  indiclal  error  oeasures  to 
evai  ate  t.ie  corresponding  reasures  for  inputs  other  than  eteps.  Relationships 
between  open-  sftd  closed-loop  forts  cf  certain  oeasures  are  also  discussed. 

Chapter  VI  discusses  -ensures  or  sensitivity  to  parameter  changes,  and 
power  rcquire^e..is,  ana  presents  the  conclusions  reached  fron  the  present  study 
c*.  criteria  for  deterministic  inputs. 
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CHAPTER  i 

CAICUIATI05  0?  DTSAHEC  P2BCOWKE 

The  dyrrclc  performance  of  &  coutrol  system  or  element  1*  assessed  by 
considering  the  following  fectors: 

1 .  Stability 

2.  Response  to  desired  Inputs 
3*  Response  to  unwanted  inputs 
&.  Accuracy 

5-  Insensitivity  to  parameter  changes 
6.  power  snd/or  energy  demands. 

The  quantitative  specification  of  dynamic  performance  consist*  of  choosing 
ne  as  urea  of  toe  above  qualities  (either  singly  or  in  caablnatioo) ,  and  setting 
desirable  values  or  linlta  upon  these  measures.  Defining  a  "peiforaance  meas¬ 
ure"  as  a  quantity  characterizing  dynamic  performance,  the  tern  •  performance 
criterion"  any  be  defined  as  a  standard  or  reference  value  of  a  performance 
measure  which  provides  a  basis  for  a  rule  or  test  by  which  ace*  aspect  of 
cynacuc  performance  is  evaluated  in  forming  a  judgment  of  »j*Un  quality. 

Per  a  performance  criterion  to  be  of  use,  it  cost  be  valid,  selective, 
and  readily  applicable.  Validity  implies  that  the  criterion  is  associated 
with  desirable  perfornsr.ee  ch&racterictico  for  the  input  environment  of  interest. 
The  requirement  of  selectivity  demands  that  the  criterion  chould  differentiate 
sharply  between  "good"  ays  teas  and  those  which  are  merely  "acceptable."  For 
the  criterion  to  be  readily  applicable,  its  expression  it.  terns  of  ays  ten 
parameters  should  be  c aspect,  and  convenient  procedures  for  its  evaluation 
should  exist. 

In  principle,  the  process  of  designing  a  rystem  to  meet  the  specified 
performance  consists  of  calculating  the  dynamic  performance,  applying  the 
measure  to  the  system  under  consideration,  and  then,  if  necessary,  modifying 
the  system  so  that  the  specified  perxermsnee  is  attained  or  approached  as 
closely  as  practicable.  For  cuuveuieucv  in  implementing  this  step-by-step 
sequence,  analytic  procedures  for  the  calculation  of  the  performance  me  a: 
should  be  available.  However,  eary  flight  control  systems  are  so  ccmplic*.  «i 
that  their  response  to  any  given  input  can  only  be  described  by  a  differential 
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equation  of  high  order,  or  by  a  large  number  of  simultaneous  differential 
equations  of  lover  order.  Physical  realities  thus  tend  to  he  ob*eur»»d  in  a 
fog  of  mathematics.  To  avoid  this  situation,  a  number  of  simplifications  can 
be  introduced  i»wo  the  analytical  representations  of  act;.al  flight  control 
systems.  Foremost  s*ong  these  is  the  assumption  of  small  perturbations  erel, 
hence,  llneerity.  The  further  assumption  is  then  commonly  made  that  changes 
in  vehicle  configuration  and  environment  occurring  during  the  notion  are 
emails  so  that  the  coefficients  of  tin.  differential  equations  are  effectively 
constant.  The  resulting  linear  constant- coefficient  equations  arc  still  of 
liigh  order.  Por  purposes  of  calculation,  a  lover-order  system  which  possesses 
(for  a  given  input)  approximately  equivalent  dominant  node  dynamics  can  be  con¬ 
structed.  This  artifice  is  particularly  valuable  in  the  calculation  of  perform¬ 
ance  criteria.  It  will  be  demonstrated  belov  that  this  core  tractable  equiva¬ 
lent  systeu  J*  relatively  easy  to  deduct  from  the  open-loop  transfer  function  of 
the  particular  loop  that  is  being  studied.  later  In  this  report  it  will  also  be 
shown  that  the  result  provides  a  sufficiently  accurate  basis  for  the  calculation 
of  performance  criteria. 

The  general  procedure  by  which  equivalent  systems  are  derived  is  *o«t 
clearly  illustrated  by  an  example.  Consider  a  pitch  control  system  for  the 
fighter  airplane  detailed  Jn  Appendix  C.  The  opcn-lcop  transfer  function  for 
the  pitch  loop  is 


Controller  Transfer 
Function 


The-  Bode  diagrao  for  G(Jcd)  ia  snovn  in  Fig.  1.  The  closed-loop  systen  haa 
three  region*  of  interest  defined  by 

|0(^)|  »  1,  over  Vhlch  |r£^|  i  1 

|<3( Jco) J  «  1,  over  vhich  1 1  |  “  end 

|g(jcd)|  :  1 

The  fona  of  the  closed  loop  transfer  function,  p~V^o^)j*  1"  this  last 
region  defines  the  "dominant  nodes"  of  the  closed -loop  systen  dynamic  response 
for  impulse  and  step  inputs.  In  cost  case6  G(J»V[l  *  r»(-Ja»)3  in  the  region 
where  Jg(<J^)|  Is  of  the  order  of  unity  can  be  a^proxicated  by  a  first-,  second-, 
or  o bird -order  systen,  the  codes  of  which  will  determine  the  cajor  features  of 
the  response.  The  open-loop  amplitude  asymptotes  cf  an  appropriate  equivalent 
ayater.  for  this  exsapJe  are  shown  in  Pig.  . . 

Applying  this  approximation  tc  the  present  example  yields  the  closed-loop 
(>»)  Bode  diagram  of  Pig.  ?.  The  Bods  ditj^rea  for  the  exact  closed-loop  system 
is  also  shown  for  cooperative  purposes.  It  will  be  observed  that  the  error  of 
th*  approximation  is  small.  If  greater  accuracy  ia  required,  more  complicated 
open-loop  equivalent  systems  can  be  produced  by  retaining  core  of  the  teres  in 
the  complete  opsn-loop  transfer  function. 

In  the  example  cited,  the-  servo  break  frequency  Is  of  the  order  of  <su.  Hare 
typically,  this  rre^yency  will  **  »  s^;  the  effect  of  the  associated  high-fle- 
quency  leads  a  nr.  lags  can  then  be  approximated  by  replacing  then  in  either  the 
open-  or  closed -locp  transfer  functions  by  a  pure  time  delay  term,  e-T3.  A  satis¬ 
factory  approximation  for  the  time  delay  is  i  «  *(?leads  *  ^lag-j)high  frequency. 
(Alternative  approximations  are  discussed  in  Chapter  IV.) 

In  general,  airplane  transfer  function  break  frequencies  and  tlae  constants 
arc  spaced  so  that  It  the  rvgion  of  crossover  can  be  satisfactorily  approxi¬ 

mated  by  a  systeo  of  not  core  than  fourth  order. 

The  equivalent  system  concept  can  be  applied  to  fora  numerical  censures  of 
each  of  the  aspects  of  dyneaic  perforsar.ee  listed  at  the  beginning  of  this 
chapter.  However,  determination  of  stability  la  usually  only  slightly  core 


—  Eounratent  system  osympto/i 
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complicated  fcr  the  cooplecc  linearized  system  coiel  then  for  the  simplified 
equivalent-  7*blc  I  cuscarizeB  techniques  far  investigating  rtability.  Since 
the  meaningful  application  of  alacxsl  all  standard  performance  criteria  requires 
tnat  ti*e  system  shall  oe  stable,  the  use  of  one  of  the  tests  listed  In  Table  I 
Cores  a  prerequisite  to  the  application  of  more  refined  criteria- 

Of  the  lemaining  aspects  of  dynamic  performance  listed  at  the  beginning 
of  this  chapter,  this  report  Is  principally  concerned  with  Item  2,  "response 
to  desired  Inputs,"  ind  the  main  body  of  the  wort  presented  here  la  concerned 
with  the  inverpretatiou  or  this  quality  in  numerical  terms  as  a  performance 
censure-  Investigation  of  the  system's  capacity  to  suppress  most  unwanted 
inputs  requires  consideration  of  gusts,  noise.  etc.,  which  are  best  described 
In  statistical  terra,  and  thus  fall  outside  the  prec*nt  report. 

However,  tcce  unwanted  incuts  (e-g.,  engine  fall’tre)  ere  veil -described  by 
deterministic  expressions;  for  these  input*  the  methods  of  the  present  report 
are  quite  applicable.  Ine  topics  of  accuracy,  insensitivity  to  parameter 
ctiangefc,  and  power/ energy  demands  are  discussed  in  Chapter  YI  In  terms  of 
perfoirance  measures- 

tonventlonally,  the  quantitative  Interpretation  of  the  response  to  desired 
inputs  is  achieved  by  scar-s  of  measures  of  the  motion  following  s.  limited 
selection  of  deterministic  inputs.  These  inputs  are: 

1.  Impulse  (Dirac  delta  function) 

2.  Step 

3.  Raa? 

.  Pover  aeries 

5-  Sine  wave. 

Besides  their  conventional  nature  as  test  inputs,  the  first  four  of  these  are 
also  representative  of  a  vile  variety  of  fll&it  control  system  '■oca and,  dis¬ 
turbance,  or  Initial  condition  inputs- 

The  well-known  convolution  relationships  of  linear  theory  {-’-6-,  Duhsael's 
integral),  spe^ieiited  techniques  described  In  fief.  3 6  and  72,  enable  Ibc 
response  to  any  deterministic  Input  to  be  calculated  fren  a  knowledge  of  the 
response  to  any  of  the  inputs  listed  above.  A  comprehensive  surrey  of  perform¬ 
ance  measures  that  have  been  prepoeed  for  deterministic  inputs  1*  presented  in 
the  following  two  chapters. 
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SAILS  1 

StttlLITY  MEASURES  *Cfl  COSSlART-COEmCTF 


SYSTEM  CWWITIES  CTYOLVSD 


direct  wasuras 

(ASSOCIATED  »&\SW«S) 


Descartes'  Rale  of  Signs 


CUractcrUtlc  equation*  o(t)  •  0 


Shw&er  of  varlatlor*  In  .ign 


rfaklal’s  Cosiitio" 


CTaracUritU:  equation*  &(e*  »  O 


Relative  value*  of  &fs)  coefficient* 


CencrallzeH  Descartes’  Rule 
of  Signs 


Modified  cMracterlatlc  equation. 

a(U  .  o  J  a(i)  - 


Huefcer  of  variations  in  tigs 


Routh-Hurvltz  Criterion 


22.6A 


Characteristic  equation,  A{s)  -  0 


Routta  test  functions; 
Kurvltz  determinants 


Cei^rallted  Aouth-Kurvltz 
Criterion 


Moll  fled  characteristic  equctlon* 

a(l)  -  0  :  £{»  . 


Ruutfc  test  functions  or  Rurvitt 
determinants  for  modified 
characteristic  equation  in  A. 


Ll"fOu--orf  Thcoreo 


Matrix  of  coefficients,  (*];  idea 
systca  equations  are  In  the  form 

[*]•  ww 


Positive  lAf*i.ltc  f&tm  or  tt<e 
jymetrlc  cat  r lx  [p] 


Mlthallov  Criterion,  Also 
[Leonhard  or  Cre^r-Leonharl 
CrlLrlso]  J5.W.:0.-^ 


Characteristic  equation  In  fore 
A( Ju>)  •  A{oj2)  ♦  joS(or) 


Koota  of  a(o£)  and  t(v£) 


Rusher  of  encirclements,  R,  of  -1.  am 
number  of  leroa,  p,  of  Kt)  -  0  In 
right  half  plane  _ 


Cauchy-Syqul st 


Opea-locp  transfer  function, 

«•>  •  4^ 


Gain  Kargin, 

* neutral  staMllty-^aetual 
Phase  cargln. 

^5*^)  *  Neutral  stability 


Peat  magnification  ratio,  Kp 


Cenc rallied  Ryquist 
Criterion  '.*9,76 


Oper-loop  transfer  function  c(s) 

.  -  -»!JO  B  1  ■  |-|  1  jV>  *  A" 


Rwber  of  encirclements  of  -1,  and 
nvesber  of  seroa  in  0(s>  vitb  real  parts 
greater  than  c,  or  damping  ratios  less 
than  (. 


Cere  rallied  gain  tmrgla, 

llStoH  -  ’ 

_ 


Oca*  rail  zed  ptese  Margin 


Precise  Root  Location 


Polynomial  Factoring 
Technique* 

27.26 


Actual  root  values  of  closed-loop 
systca:  [j,  oj,  l/?t. 


2.  Root  Locus 


J.  Generalized  C(s) 


Characteristic  equation,  A(s)  •  0 


Poles,  zeros,  and  gains  of  open-loop 
tr«*afcr  function.  C(u) 


Open-loop  transfer  functions  G(s) 
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TABUS  Z 

.^ABILITf  HASdRES  TOR  COSSTAW-COWIOIZST  LINEAR  bfSTQe 


APPLICATION  nrawo^s  ASD  AID6;  RZMAftth 


nr.sc?  JStsuRir 

ASSOCIATED  ffiASUTES) 


P.u'.th  test  Fun Tticui  or  Kurvltz 
Actcmina.it*  for  nidified 
:.*4?ict«rlsiic  cq-ntloa  In  X. 


Positive  definite  fora  o'  toe 
eyrcetric  mtrix  w 


hoots  of  A(u*0  and  B(ur) 


All  coefficient*  of  ssae  algn  «re  o  cectmry,  but  not  sufficient,  condition  for 
stability.  <iwwwi  of  yoiltirt  £ne«a*lv«]  real  root*  cannot  exceed  tbe  maker  of 
variations  of  sign  of  A(a),  [&(-s)j. 


All  coef-i'l'ut*  of  aare  sign  are  a  necessary,  out  uot  auffldcnt,  condition  for  all 
roots  to  have  negative  real  parts  Ices  than  assigned  v*li«  of  -o. 


Provide  necessary  and  sufficient  conditions  for  stability;  also  determines  number  of 
roots  in  right  half  plane.  Routb'c  Algorithm  Is  a  staple  aid  for  development  of 
hlcher-c-nter  test  fu.*ctlcas.  One  test  function  is  critical  when  a  paraneter  variation 
causes  a  cnange  from  stability  to  Instability.  (Ref.  22) 


Provide  necessary  and  sufficient  conditions  for  roots  to  have  real  part*  less  than 
assigned  values  of  -c- 


S  vs  ten  is  asymptotical!*  stable  If  tbe  matrix  [PJ,  which  Is  svnetrlc 
([p]  •  Transpose  [PJ  •  tP1]  ),  satisfying  tbe  aatrlx  equation  [aIM  ♦  [p}[a]  »  -[l] 
(toe  unit  ratrlx)  is  positive  definite  (all  principal  minors  cf^PJar*  positive),  the 
cl'rentr  of  £?]  are  found  frea  the  nfn  ♦  \)/2  elmultereoua  equations  resulting  free 
expending  the  tmtrlx  equation  above.  The  theoren  can  also  be  ateted  in  tens  of  the 
Uapounoff  fwictlon  Y((:])  * 


Type  and  relative  sequence  of  roots  indicates  condition  of  stability  or  instability. 
In  stable  situation  roots  of  A(c£)  -  0  and  B(u£)  -  0  are  staple,  real,  and  positive; 
and  they  alternate  in  tb~  sequence:  ABA  B,  etc.;  usually  applied  graphically. 


mmmsmm 

ClTca  number  of  zeros,  Z  «  P  ♦  a,  of  1  •  C(s)  -  0  in  right  half  plane.  Ordinarily 
applied  graphically  v".-  polar  plot  or  G(wVt)  Bode  dlagraas. 

Geir.  r'-gln, 

neutral  slot ility^actval 

Gain  change  necessary  for  rza’rvl  C* ability. 

Phase  margin, 

9(afc)  -  Neutral  stability 

r««  vMH<ec  n  .paired  for  neutral  stability,  with  gain  held  constant. 

Pea*  sogni*i cation  ratio,  Kp 

Measures  =axl**»  eicsed-loop  resonance.  Usually  determined  fro-  open-loop  plots  (polar 
or  logarlthaie  gain-phase)  and  closed-loop  nvarlaya  (M  circles  or  Rlcbols  chart). 

I 

Busber  of  encirclement*  of  -1,  and 
oust*  of  zeros  in  5(s)  with  real  pets 
greater  tnan  a,  or  dacplAg  ratio*  less 
than  1 . 

Circa  nvsber  of  zero*  of  1  «  c(s)  •  0  which  have  real  parts  greater  than  r,  or  damping 
ratios  less  than  {•  Ordinarily  applied  with  C(a)  Bods  diagrams. 

I 

Generalized  gain  aorgin, 

U°W1  • 1 

''actual 

Calr  change  required  to  achieve  roots  with  specified  o  er(. 

■ 

Generalized  phase  aargln 

™ "*■'  '■•t-f.?-nu!.wwvniani.s.'.'gi  nimu  u  ».-■,<  vmi  nmuM 

■ 

Actual  root  values  of  eloscd-loop 
systex:  Ct,  a*,  ,/?i- 

Coeplete  def initio.-.  oT  system  transfer  characteristics,  with  knowledge  of  Input, 
response  Is  completely  defined. 

1 

Techniques  Include:  Ncvton's  method,  Berner's  aetbod,  synthetic  division,  Graeffe's 
root-squaring  method,  Ida's  method,  etc. 

1 

Supplementary  Techniques  In  Unified  Servo  X-talysl*.  (Ref.  JJ) 

2 

CHAPT3F.  rr 
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The  tern  “ay t tea  characteristics"  will  be  used  to  denote  all  performance 
seasures  not  expressed  as  explicit  functions  of  error.  For  example,  phase  car- 
gin  *M  ttcs-io-peck  arc  both  ays  tea  characteristics,  whereas  ITAE  is  classi¬ 
fied  ad  en  "i**Iicial  error  aeasure"  and  is  discussed  in  Chapter  HI*  Although 
an  infinite  nueber  of  systes  characteristics  could  be  devised,  useful  performance 
erasures  in  practice  are  obtained  only  fron  characteristics  directly  descrijn.  /e 
of  the  ays ten  transfer  function  or  the  impulsive  or  indicial  (step-input)  re¬ 
sponse. 

Tables  H-A  and  II-B  susxsarlre  a  broad  cross-section  of  systes  character¬ 
istics.  All  the  characteristics  that  have  been  found  in  the  references  listed 
at  the  end  of  tMs  report  are  included,  wish  one  exception  (the  product  of  peak 
overshoot  and  tine -to- peak),  which  is  discussed  on  page  27-  The  characteristics 
listed  in  Table  H-A  are  quantities  directly  obtainable  froe  G(s),  G(jco),  or  the 
closed  loop  fores  1  It  is  possible  to  c^r-emllse  sore  of  the 

G(jo.)  measures,  such  as  phase  rargin  and  gain  margin,  into  analogous  quantities 
for  G(s)  Terrs.  Apart  free  the  error  coefficients,  the  remaining  entries  on 
Table  n-B  describe  the  response  to  a  step  input.,  and  will  be  called  "iniicial 
response  characteristics." 

Ideally,  the  irxlicial  response  characteristics  wild  partition  the  response 
into  the  regions  indicated  in  Fig.  >.  "Dead  tire"  in  the  time  to  attain  10  per¬ 
cent  of  the  final  value,  rise  time  is  the  tire  to  go  free  10  percent  to  90  per¬ 
cent,  and  decay  tire  is  the  time  for  the  transient  to  decay  free  90  percent  to 
within  3  percent  of  the  final  value.  The  sun  of  the  dead  tire,  rise  tire,  and 
decay  tire  la  known  as  the  settling  tire  (which  can  be  defined  directly  as  the 
tire  for  the  indicia!  response  to  reach  and  thereafter  regain  within  3  percent 
of  the  final  value) .  Unfortunately,  only  the  settling  and  rise  tires  are  read¬ 
ily  related  to  transfer  function  characteristics,  so  the  quantities  involved  in 
the  Idealised  partitioning  of  the  indicia)  response  are  replaced  by  the  soce- 
what  overlapping  reasures  indicated  on  Table  1I-B- 
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In  order  to  assess  the  ce.*it  of  any  given  system  dynanlc  performance  from 
its  indicial  response,  it  is  necessary  to  define  "good"  dynamic  performance  in 
terms  of  the  indicial  response.  It  is  generally  agreed  that  "good"  dynamic 
performs implies  lov  overshoot,  short  dead  tine,  fast  rise  time,  and  good 
damping  of  motions  subsequent  to  the  decay  time.  (The  last  requirement  implies 
lov  settling  tine.)  Thus  a  good  indicial  response  will  comply  with  certain 
specifications  on  its  "shape,"  which  can  be  defined  in  terms  of  overshoot, 
ratios  of  dead-time  to  settling-time  and  rise-time,  etc.  However  specifica¬ 
tions  of  "shape"  alone  is  insufficient  to  ensure  that  the  indicial  response 
will  be  satisfactory;  some  parameter  defining  the  time  scale  of  the  response 
rust  also  be  specified.  The  Implications  of  this  last  requirement  will  now  be 
discussM  briefly. 

The  over-all  time  scale  of  the  indicial  response  depends  upon  the  bandwidth 
of  the  system.  Practical  considerations  of  inertia,  power  demands,  etc.,  result 
in  increasing  penalties  in  weight  and  complication  as  bandwidth  is  increased. 
However,  In  this  genualized  investigation  it  is  not  possible  to  set  these  upper 
limits  upon  bandwidth  explicitly.  For  statistically  described  Inputs  and  for 
deterministic  inputs  (such  as  rectangular  pulses)  which  arc  of  finite  specified 
durations,  lower  limits  on  bandwidth  can  be  set  at  least  approximately.  For 
example,  the  settling  time  should  not  exceed  the  pulse  duration.  However,  for 
impulses  end  step  inputs,  upper  limits  on  settling  time  cannot  te  specified  in 
the  absence  of  further  information  regarding  the  operating  envlrowaent.  Thus 
this  chapter  and  Chapter  HI  are  essentially  limited  to  a  study  of  those  aspects 
cf  dynamic  performance  which  can  be  represented  by  the  "shape"  of  the  Indicial 
response.  To  focus  attention  oa  "shape*  rather  than  time  scale,  performance 
measures  such  as  settling-time,  rise- time,  etc.,  are  all  expressed  in  nondiaen- 
sional  forms.  For  example,  second-order  system  characteristics  such  as  rise- 
time,  settling-time,  etc.,  arc  normalized  through  multiplication  by  a*),  where 
is  the  system  undamped  natural  frequency.  A  more  general  procedure  for  nor¬ 
malization  will  b.-»  given  at  the  beginning  of  Chapter  HI. 

The  general  procedure  for  obtaining  system  cliaractcristics  is  to  construct 
the  appropriate  transfer  function  representation  or  indicial  response,  and  read 
the  measure  directly.  Thus,  crossover  frequency,  phase  and  gain  margins,  and 
frequency  of  instability  nay  be  obtained  free  the  open-loop  C(.Vo)  Bode  diagram. 
Application  of  v^e  unified  servo  analysis  method  of  Ref.  55  (hereafter  referred 
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to  as  U.S.A.M.)  facilitates  swift  construction  of  the  closed-locp  G(>c)  Bode 
diagram:,  inspection  of  which  yields  peak  magnification  ratio  and  frequency,  and 
bandwidth. 


For  second-order  ays  tecs,  it  is  poeaible  to  develop  <*xact  formulae  for 
cany  system  characteristics  and  approximations  for  the  remainder.  Tables  II -A 
and  II-B  list  these  formulae  in  terns  of  open-  and  doeed-loop  parameters  for 
a  system  having  the  open-loop  transfer  function. 
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and  a  closed-loop  transfer  function, 
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whirrs  is  the  unda-ped  natural  frequency,  and  £  is  the  daaplng  ratio.  The 
closed-  and  open-loop  parameters  are  related  by  the  equations 
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(This  second-order  system  corresponds  to  a  wide  variety  of  equivalent  systems 
encountered  in  flight  control  applications) .  Elementary  operations  on  these 
relations  enable  exact  formulae  to  be  found  for  all  transfer  function  frequency 
domain  characteristics.  The  resulting  formulae  are  listed  In  Table  H-A,  to¬ 
gether  with  references  where  derivations  can  be  found.  Simpler  approximate 
formulae  say  sometimes  be  preferred,  and  two  standard  approximations  are  dis¬ 
cussed  below. 


Figure  4  Illustrates  a  linear  approximation  to  phase  nargln 
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Cccpericon  with  the  exact  result  indicates  acceptable  accuracy  for  <  .*>0  deg. 
Bandvidth  nay  also  be  found  approxlcatcly  using 


CV  '  0  ♦ 


(7) 


which  is  obtained  by  neglecting  terse  in  £  in  the  exact  expression  in  Table  XU. 
Equation  7  io  coopered  with  the  exact  band  vines  fer  a  second-order  system  in 
Fig.  5- 

Pbaso  uargin,  bandwidth,  and  peak  signification  ratio  are  widely  used  for 
performance  apecifi cation.  The  phase  margins  of  all  the  standard  forms  pre- 
eented  in  Chapter  HI  (Table  VH)  are  between  50  and  (0  deg.  In  systems  with 
a  detain  ant  second-order  node,  this  would  be  expected  to  yield  ndc  quite  (£  i  0*7) 
clasping.  Similarly,  elird. nation  of  the  frequency  response  peak  leads  to  ade¬ 
quate  damping  of  the  dominant  nodes-  The  connection  of  these  measures  with  the 
transient  response  is,  in  general,  neither  unique  nor  explicit,  except  for 
occur.!  eric.  „„  .. .  discvn»»eu  iu  caaaeclion  with  Table  111  - 

A  good  system  for  a  given  application  must  or  necessity  hare  phase  margin, 
bandwidth,  and  peak  magnification  values  which  lie  within  relatively  narrow 
limits,  but  a  system  which  complies  with  these  limits  is  not  necessarily  good- 
This  follows,  or  course,  from  the  fact  tnat  the  behavior  of  the  actual  transfer 
functions  of  Interest  is  defined  only  in  a  gross  ccncc  over  a  narrow  frequency 
band  by  these  particular  measures.  In  general,  therefore,  none  of  these 
characteristics  taken  alone  yields  a  valid,  selective,  and  reliable  measure. 

Of  the  indici&l  response  characteristics  listed  in  Table  U-3,  exact 
formulae  exist  only  for  tlae-to-pea*  and  peak  overshoot  of  second-order  systems. 
For  higher-order  systems,  time -to- peak,  settling  time,  and  rise  time  can  be 
estimated  for  so me  special  classes  of  systems  by  means  of  the  charts  of  Ref.  13, 
17,  and  2^.  Reference  is  presents  rise  times  for  several  classes  of  third- 
order  systems;  the  results  are  discussed  on  page  Jl.  The  rest  of  this  chapter 
is  mainly  devoted  to  a  discussion  of  the  calculation  and  interrelation  of 
second-order  system  characteristics. 
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TABLE  HI 


ATION  OF  SECOND-ORDER  UNIT-NUMERATOR  SYSTEM  CHARACTERISTICS 
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FORMUIAE  FOR  IHIIIC2AL  RESPOuSE  OiARACinttSTICS 


Exact  fcn^luc  arc  not  available  for  scoe  of  the  indlclal  response  charac¬ 
teristics,  and  the  range  cf  validity  of  the  approximations  quoted  in  Table  II -B 
requires  some  consideration. 


The  indiclal  response  of  a  closed-loop  system  haring  the  open-loop  trans¬ 
fer  function  c(s)  »  ^  ^  xs 


where 


c(t)  -  1  ♦  /  '  e'ten1  «in  (o^  >/l  -  C2  t  -  V) 

V'  -  r 


-C 


t2 

— »  <  -  cos  C 


(8) 
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This  equation  can  be  solved  eraetly  for  •  0,  yielding  tine-to-peek  and 

peak  overs-oo  . 

Settling  tine  is  cccmcnly  approximated  hy  considering  only  the  envelope 
of  the  response 
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(5) 


For  an  overshooting  response  c(t#)  -  l.0>.  Putting  ce(t$)  »  1 .05  in  Eq  9 
yields  . 

.  do 

Equation  10  differs  free  the  approximation  given  on  pages  22-hl  of  Ref.  36  where 
the  -  C2  bern  is  replaced  by  unity.  Figure  6  compares  Eq  10  and  tnc  further 
approximation  ts  •  j|-  given  in  Ref.  1  with  the  settling  time  obtained  by  direct 
measurement  of  analog^ computer  responses.  It  shoild  be  noted  that  approximations 
to  settling  time  (e-g.,  do  not  reproduce  the  sawtooth  shape  of  the  exact 
graph. 
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The  Equivalent  TL-*:  Constant  measure  is  m03t  appropriate  for  systems  .'ith 
a  'Iccirant  first-order  node.  To  toe  extent  toat  this  obtains,  Teq  i  1  A*.  asd 
n  particularly  simple  connection  exists  bel*c<.n  the  frequency  o-.i  time  dcrains. 
When  extended  to  second-order  systems,  the  variation  of  Equivalent  Tim*.  Constant 
with  C  in  more  regular  than  the  settling  time,  although  the  measure  is  actually 
inappropriate.  For  instance,  Fig.  7  i  Dust  rates  measured  ?eq,  and  the  approxi¬ 
mation  of  Ref.  >5  »  i/«^  for  a  at-cond-order  system,  with  the  open-loop  trans¬ 

fer  function  of  Eq  2.  Rote  that  the  minimum  value  of  Tcq  occurs  at  \  *»  0. 

Delay  tine  is  defined  here  as  the  tine  for  the  indiclal  response  to  % 
achieve  £0  percent  of  its  final  value.  The  exact  value  is  ccmpared  vith  tre 
appro-'i-ation  (given  ir.  Ref.  72) 


In  Fig-  3.  This  approximation  is  satisfactory  only  in  toe  "opt Inal"  region  of 
0.6  <  i<  O.y.  An  empirical  linear  relation  tj  ■  is  more  generally 

applicable.  Again,  as  with  Te~,  minim.  I-.  delay  tine  is  achieved  at  £  »  0- 
Equation  II  is  extents!  to  higner-order  systems  in  Chapter  IV,  where  it  is 
shown  that  for  optimal  systems,  simple  and  accurate  approximations  to  delay 
tine  can  be  obtained. 


Rise  tine  is  defined  as  the  time  for  the  indicia!  response  to  rise  Tree 
10  to  rO  percent  of  its  final  value.  Reference  *6  presents  the  simple  aj> 

j  a 

proximation  •  This  is  compared  with  the  exact  rise  tine,  and  a  more 

refined  approximation  (Sq  12)  in  Fig.  - 
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A  minimum  rise  tine  system  possesses  a  lo-  £.  An  approximate  formula  for 
rise  time  is  given  is  Ref.  25 
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where  Ej  ard  sre  the  velocity  and  acceleration  dynamic  error  coefficients 
defined  at  the  bottom  of  Tabic  U-3.  For  the  second  order  syste-  considered, 
the  accuracy  of  this  approximation  is  far  inferior  to  Eq  12. 


25 


[ 


Figureu  10  and  11  Illustrate  IN'  7“sk  o«*“r*hoot  and  time-to-peak,  respec¬ 
tively,  for  a  unit -numerator  second-order  system.  The  minima  occur  at  £  -  i  and 
?-0,  respectixely. 

It  should  now  be  apparent  that  a  criterion  based  on  a  single  indicial  re¬ 
sponse  time  or  transfer  function  measure  is  unlikely  to  be  valid  for  the  variety 
of  systems  encountered  in  flight  control  optimization.  Most  of  the  measures 
specify  only  particular  elements  of  the  response,  *jri /or  are  restricted  to  spe¬ 
cial  tyres  of  G(s)  behavior  in  the  region  |g(>o)|  *  ’ .  Minimization  of  the 
measures  results  cither  in  low  damping  ratioo,  io«  i-Jume  "yrglns,  etc.,  o**  is 
unselective.  The  or  u3e  of  the  indicial  response  and  elementary  transfer 
function  measures  is  for  specification  purposes;  this  will  be  discussed  further 
bclc*-. 

One  indicial  response  measure,  i.e.,  settling  time,  is  a  possible  exception 
to  the  general  state-sent  above.  Minimum,  or  nearly  minimum  settling  time  is 
often  used  no  a  criterion,  and  because  it  gives  a  £  »  0.69  for  second-order  sys¬ 
tems,  it  is  worth  examining  is  the  light  o*  r**jiir*»*nt«s  tor  '*ri*'*»*ia. 

Settling  tine  is  a  possible  exception  to  tl.e  tNtateae.it  above.  Since  "good" 
indicial  response  implies  low  settling  time,  it  is  reasonable  to  inquire  whether 
tne  eiiUrisn  of  low  settling  time  always  results  in  indicial  responses  which 
are  satisfactory  in  relation  to  overshoot  and  the  other  parameters  defining  the 
"shape*'  of  the  response.  Figure  6  shows  that  for  a  unit-numerator  normalized 
second-order  system,  tne  criterion  or  billing  time  selects  £  «  O.69. 

However,  the  sawtooth  graph  yields  almost  equal  settling  times  for  0.*»3  <  £  <  0.60, 
with  sudden  Jumps  at  each  end  of  this  range.  Similar  discontinuous  selectivity 
characteristics  occur  with  highei -order  syst.*as,  as  shown  in  Ref.  32  (which  pre¬ 
sents  an  exhaustive  investigation  of  settling  time 3  for  higher-order  systems, 
including  standard  forms  of  first-  through  eighth-erder) •  These  discontinuities 
can  be  removed  by  using  approximations  to  settling  time  (such  as  tft  »  3/£o*»  see 
Fig.  6),  but  ruch  approximations  re  difficult  to  derive  for  higher- order  eye  terns, 
uiuc. 5  dominant  modes  can  be  identified.  In  general,  therefore,  minimum  settling 
tine  cannot  be  reccraaenaed  for  use  as  a  sole  criterion. 

Minimization  of  the  product  of  time-to-peak  and  peak  overshoot,  proposed  as 
fc  criterion  in  Ref.  £*,  is  also  examined  in  Ref.  32.  Its  failure  to  discriminate 
against  undershooting  responses  and  frequent  selection  of  very  poor  response  moke 
this  characteristic  unsuitable  for  use  as  a  performance  criterion. 
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The  oc c  of  error  inefficient*  performance  criteria  is  discussed  in 
Chapter  V-  It  is  shown  that  each  error  coefficient  ir  pi^porticcail  to  a  l;ac- 
weighted  integral  of  the  error  response  to  an  lspulei**e  injait.  The  error  coef¬ 
ficients  are  thus  the  sus  of  positive  tod  negative  error  integrals,  and  fail  to 
indicate  whether  a  anal!  value  of  this  sum  represents  a  snail  absolute  error, 
or  a  small  difference  between  large  positive  and  negative  error} *  It  is  con¬ 
cluded,  therefore,  that  err«*  ,'~*r  'icientc  by  themselves  do  not  provide  valid 
performance  criteria  for  general  systems. 

GENERAL  CCRRBIATI03  07  SYSTEM  CKARACTnOOTICS 

It  is  sot  surprising  that  none  of  the  indioial  response  characteristics 
considered  in  this  chapter  provides  an  acceptable  performance  criterion.  These 
characteristics  are  suitable  lor  performance  specification,  rather  than  opti¬ 
mization.  For  example,  it  my  be  convenient  to  specify  system  performance  in 
terms  of  bandwidth  and  tine-to-pcaX.  It  is  important  that  any  specification 
framed  in  terms  of  system  chn racte ri  *  v  i  c*  abode  be  practically  realizable  (l.e., 
mutually  exclusive  value «  or  Units  of  characteristics  must  be  avoided).  To 
facilitate  this  process,  it  is  desirable  tc  have  o  complete  correlation  express¬ 
ing  each  system  characteristic  in  teres  of  any  other  ays  tea  characteristic. 

Table  HI  bar  been  prepared  for  this  jurpoae.  Crossover  frequency,  bandwidth, 
phase  margin,  peak  frequency,  magnification  ratio,  tine-to-peak,  peak  overshoot, 
and  delay  time  arc  all  expressed  in  teres  ot  <r^,  and  each  other,  for  a  second - 
order  rye tern. 

For  flight  control  system,  despite  the  presence  of  nuaerct*  unalterable 
elccerts  eetc-iited  the  vehicle  configuration,  the  variety  of  possible 
system.  3  c  so  groat  that  no  measure  of  the  loiici&l  response  based  upon  any 
single  instant,  or  measure  of  transfer  function  characteristics  at  a  single 
frequency,  can  hope  to  provide  wore  than  necessary  rather  than  sufficient  con¬ 
ditions  for  system  goodness.  Recognition  of  this  fact  has  led  cany  investiga¬ 
tors  to  propose  assures  based  upon  integrated  functions  of  the  indicia!  error 
response.  These  Indicia!  error  res* urea  air  discussed  in  th~  next  chapter. 

The  preceding  discussion  of  system  characteristics  has  examined  the  merit 
of  each  characteristic  a*  s  performance  criterion  for  unit-nuasretor  second- 
order  systems.  Ttls  type  of  syatea  vxs  selected  because  mny  equivalent  flight 
centred,  systems  are  in  this  category,  and  consequently  it  provides  a  fair  basis 
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for  lcltld  assessment.  Systen  charecteri&tics  that  do  rx>t  yield  satisfactory 
■jcrforEar.ee  criteria  for  secood-order  sv*^.*tx  reed  not  be  considered  further  In 
the  tearch  for  a  satisfactory  performance  criterion-  The  diversity  of  possible 
third,  and  hlgher-ordci  &yt»v«*s  precludes  a  general  correlation  of  systca  char¬ 
acteristics,  as  has  been  achieved  for  secood-order  systems.  However,  taere  is 
case  value  in  collecting  the  Halted  data  available  od  syotea  characteristics  for 
third-  and  higher-order  systems,  and  ccs-paring  tnc  results,  vi.urv.vcr  practicable. 

The  well-known  charts  of  Chestnut  and  Mayer  (Ref.  13)  reproduced  In  Ref.  j6 
correlate  tiae-to-peak  and  settling  tine  with  frequency  response  characteristic*. 
Elgerd  (Ref.  2*0  presents  Icdlclal  response  tine  histories  of  a  variety  of  third* 
orW  systems  from  vnich  rise  tine,  tim-to-peaX,  etc.,  cay  be  nutw-reu  directly. 
To  examine  the  implications  of  his  results,  and  to  cospare  then  with  those  ob¬ 
tained  from  other  sources,  a  standardized  third- order  unit- numerator  system  is 


m  (■-  >'?)[**<  >•!*<  (!•««*] 

This  fera  <«  not  cov^rad  bv  the  charts  of  Ref.  36,  but  Is  discussed  by  Clenent 
(Ref.  1 7) ,  vr.o  presents  rlcc  tins,  settling  tine,  and  peek  overshoot  for  * 
norcalized  3yetem  having  the  following  transfer  function: 


^  (1  ♦  6s)(e2  ♦  2£s  *  1) 

figure  If  illustrates  Cleservt's  results  on  settling  aloe  for  tnc  standardized 
system  of  Sq  1*.  The  points  aerked  dooc^e  val»s  chewed  by  taesns  of  the  tran¬ 
sient  response i  of  Ref.  2^.  The  agreement  is  goad,  although  allowance  must  be 
cade  for  the  fact  that  the  sawtooth  shape  of  tlo  settling  time  graph  has  been 
scoot  bed  in  Ref.  17- 

Burnett  and  Shumate  (Ref.  13)  concUte  rise  tire  with  peak  power  for  a 
variety  of  thl’-d- order  systems.  TV;  resulting  values  of  rise  tiae  and  over¬ 
shoe-.  for  the  systes  of  Eq  are  presented  in  Fig.  13,  from  which  it  will  be 
seen  that  the  agreement  with  Ref.  2*  is  t;cncrally  close. 


CHAPTER  III 

IKDICIAL  ERROR  MEASURES 


The  ter*  "indlclal  error  Measure"  will  be  used  to  denote  integrated  func¬ 
tion?  of  the  error  response  to  a  step  input.  Table  IV  suamarlzes  all  the 
indlcial  error  Measures  that  have  been  suggested  as  criteria  in  the  references 
listed  at  the  end  of  this  report.  Most  indlcial  error  Measures  are  of  the  fora 


o  f00  - 
ITH!  -  j  t  *U4t 

•'o 


etc.. 


06) 

07) 


(18) 


where  U  is  a  general  function  of  error,  cuch  os  E2,  |e|,  etc.,  not  involving 
tine  explicitly. 

These  measures  are  meaningful  only  for  zero-position-error  systems;  most  of 
the  published  investigations  of  the  relative  Merits  of  these  measures  have  been 
concerned  with  their  application  to  systeas  having  unit-numerator  closed- loop 
input -output  transfer  functions.  Since  many  flight  control  systems  reduce  to 
unit-fuse r» lor  equivalent  systeas,  this  is  a  fair  as  well  as  convenient  basis 
for  assessment  of  criteria.  To  achieve  compact  presentation,  it  is  customary  to 
normalize  the  systems  and  associated  performance  seasons.  Thus,  a  cystem  having 
the  transfer  function 

.  ...  ♦  PgT2  .  P)7  ♦  Pp 

*  •••  ♦  V2  *  V  *  «c 


where  7  is  the  Laplace  transform  vari¬ 
able  at  this  point  (elsewhere  s  is 
employed) 

can  be  put  in  normalized  form  by  the  procedure  of  Ref.  37  >  reproduced  below  for 
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*\f  reference. 


1 .  Define  a  constant  Dg  so  that 


fio 


% 

«o 


(20) 


2.  Define  new  coefficients  for  the  nuaerator  and  denominator  terns 


Q< 

4l  “  »  1  -  '•  2>  "•  B 

(21) 

«n 

?1 

1«  »  — — : —  »  1  a  C,  2,  a 

1 

(22) 

3.  Divide  the  numerator  end  denominator  of  Eq  19  by  and  apply  the 
definitions  of  Eq  20,  21 ,  and  22.  The  transfer  function  then  be  come  a 

c(r)  2a°oE'V  *  —  ♦  *  PiPq*'1*  *  fppg 

^  '  r”  *  ViV”*’  ♦  •«  ♦  <&r2?  ♦  ♦  «8 

k.  Introduce  a  new  transfer  variable  so  that 


s 


JL 

»0 


(2>0 


Tr.t-.n  the  transfer  function  reduces  finally  to  the  noraalired  fora 


C(s) 

R(«) 


•»p2s  +Pt8+P0 


Vic 


(25) 


Froe  diaensioral  considerations  it  cor.  to  shovn  i4ml  to  convert  a  novrsalized  per¬ 
formance  measure  to  its  denorsalized  form,  the  following  relations  should  be  used: 


:!,<v> 

■  IU(ncr-’> 

(26) 

‘  iraW» 

(27) 

it2u(^» 

- 

(23) 

The  normalized  fora  of  the  closed-loop  sec.'ui-uraer  system  con>iu<i-mi  previously 
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TABLE  IV 

1XDICIAL  ERBGR  keasures 


— 

IOH 

_ , _ 

OOTERIOa  ASSESSMafT 
(EASED  UPON  UKIT-KUJGlATCfi  SXSCTe) 

ASSOCIATED  APPU  CATION  TBCHKiaJES  AND  UDSj  SaJAHtS 

■ 

If  conn  trained  to  ncnovershooting  second-order  sy&lca, 
t  -  1 .0.  If  overshoot*  -re  allowed,  £  -  o  is 
eeiected.  lack  cf  validity  and  selectivity  eliminates 
as  a  criterion.  (Kef.  57) 

For  impulsive  response  (not  step)  each  of  these  measures 
is  proportion!  to  a  particular  error  coefficient. 

(See  T»ole  XX-B  and  Chapter  V) 

If  constrained  to  nonovershooting  seccnd-ov»W  system, 

£  •  l .0.  Applies  nr a vie r  weighting  (ard  wore 
selectivity)  to  system  *;tii  £  >  1  ******  t:«>  control 
area  criterion.  If  overshoots  ore  allowed.  (  «  0  is 
selected.  Inc  a  of  general  validity  eliminates  as  a 
criterion.  (Kef.  37) 

■ 

.oo 

Similar  to  J2  t«(t)dt,  except  weights  responses  for 

5  >  1  even  mere. 

1 

For  norcalized  second-order  system,  criterion  selects 
£  »-  0.63.  Criterion  is  moderately  selective  on 
low -order  unit  numerator  system,  but  nonselective  on 
zero- velocity  error  systems  and  higher-order  systems. 
Easily  mechanized  on  analog  computer.  Lack  of 
selectivity  eliminates  ae  a  criterion.  (Ref.  37) 

Analytic  form  given  In  Chapter  in  lEq  55) 
and  Appendix  A. 

I 

Selective,  reliable,  and  easy  to  apply  for  all 
sys tests  (through  eighth-order)  investigated. 

Approaches  the  ideal  criterion  for  routine 
optimization  employing  analog  computers.  (Ref.  37) 

Criterion  is  thoroughly  explored  and  la  good  in  nearly 
every  respect,  except  complicated  nature  of  analytic 
fores »  Analog  computer-derived  standard  forms  and 
transient  responses  exist  for  zero-displacement  error 
systems  through  eighth-order,  zero-velocity  and  zero- 
acceleration  error  systems  through  sixth-order. 

(R«f-  57) 

■ 

Highly  selective,  giving  £  -  0.6  for  a  normalized 
second-order  unit-numerator  system.  Analytic  fora 
is  very  complicated  and  less  easy  to  apply  than  ITAE; 
also  less  convenient  than  ITAE  for  analog  computer 
mechanization. 

Criterion  is  much  leas  thoroughly  explored  than  ITAE. 
Possibly  not  worth  the  extra  complication  over  ITAE. 

Selects  a  value  of  £  ■  0.5  for  second-order  systems, 
lack  of  selectivity  and  specification  of  highly 
oscillatory  responses  for  higher-order  systems 
sakto  use  os  a  criterion  unsuitable.  (Ref.  37) 

Is  the  simplest  of  the  higher-order  measures  to  apply 
analytically  (Ref.  63),  viz: 

•C'2(t,dt  .  jjj  i{-s)E(.)d, 

which  is  a  well-tabulated  integral.  (Ref.  bb,  38) 

General  and  optimal  (standard)  fonss  and  transient 
responses  through  fourth-order  exist  (Tables  V  and  VI) 
and  allow  a  fairly  simple  detr  mi  nation  of  effects  of 
uncertainties  in  system  parameters.  (Can  be  evaluated 
fairly  easily  by  operations  on  normal  desip)  charts.) 
These  two  advantages  make  its  use  as  a  criterion  cowan 
in  Spite  of  its  relative  lack  of  selectivity  and 
validity. 

2 


J7 


/ 


TAW*  IV 


(exartiaun>) 


sweox 

HtUSUltZ 

a.'^socia'.td  carmuc* 

rrr2 

u?Mst 

Klnbrjs  value70'75 

r#1«u  »  f  •  o.’fA  tor  »»co«d-ord*r  eyaten. 
(Rtf.  i7 )  Kw  fair  selw.Ivlljr  vd 
reliability  for  hi*ner-crder  systems 
(through  fourth  order).  {Ref.  70)  Appear* 
very  promising  aa  aa  Ideal  unitary  criterion 
In  analytic  studies. 

n”E= 

TO.Ty 

Minima  value  1 

Tor  9scood»crder  ayrtaas,  5  1  0.65  for 
n  •  2-  (Raf.  J?)  Selectivity  become 
greater  as  n  Inert sees,  and  criterion  also 
becoars  sore  complex  to  assess. 

5 

TIHds  5  •  0.A67  for  seccnd«order  system 
for  «—»  co.  (Ref.  73)  Kss  not  yst  teen 
evaluated  for  Mgnsr  ord*’>  eystecs. 

/•li-pwi* 

WhU*  another  pareocter,  ilka 

J ^  (tej-que/^dt  «  const. 

Tails  the  resdy  applicability  test  for  the 
constant  coefficient  systems  considered 
here.  Choice  of  »  ti  sxvltrsry,  asking  the 
general  fora  of  this  Measure  difficult  to 
assess  In  the  present  t«r-s  of  reference. 

>»»  v««  60-e' 

ProUilj  falls  ready  applicability  test. 
Ounce  or  >  is  erottrery,  sssing  ta#  gsner&l 
font  of  this  Mature  111 -define!  and 
difficult  to  arses*. 

1 

J? ['<*-)  -  -(«■>•■<*- !„)f « 

vtwre  i»/t  la  *  wit  *»*j* 
erd  T©  la  *  tlaa  deity 

KllU  value  ^ 

ASSOCIATED  HSPORSS  Wa/atS 

Sot  yet  evaluated.  A  logical  font  for 
syst-ea*  with  finite  position  error.  Choice 
of  T0  1*  arbitrary. 

1 

-  *0>  -  *(*>!*« 

.  ,  66,09 

Klalcua  value 

Sot  yet  evaluate*,  lack  cf  fixed  value  fer 
?0  aakee  this  eeesure  vague  and  laic finite. 

(ikjmW^  Input) 

tac{t)dt  (Upulae  laput) 

« -X”  c<t>4t  5 

s  -X° «MtJ» 

vbers  »  iif  predetermined 
cmslttt 

Unactisractcry  response  for  higher  order 
systems.  (Ref.  3*) 

1  * 


TAALfc  IV 

(ccrnxuEJ) 


ASSXifcru  carrouca 

WITMCa  A SSE3SH0G 

{bash)  urcs  uxrr-xMJATC»  srsrnc) 

AXXlRTTf  AimCAtlOK  mwifcfcS  ATS  A7I~,  IlCWOX 8 

,  111,7/ 

'•cud  v*au« 

Select*  «  {  *  0.59*  for  second-order  tyitn. 
(Ref  yj )  Hus  fair  aele-  vity  vA 
reliability  for  hlgher-crder  ijnUsi 
(through  fourth  order).  (Ref.  70)  Appear* 
very  proofing  *t  an  Ideal  unitary  criterion 
.n  analytic  studies. 

Relatively  easy  to  apply  analytically  as  it  bo  evaluated 

using  reae  Integral  fora  as  noted  above,  vie: 

L  -  ■  .if  0  3;  jej  /,.X*!s'')I(e'’)45} 

vhere  0  la  carried  as  a  paramter  in  the  Integration.  General 
and  optlaal  (standard)  fore*  and  transient  responses  through 
fourth  ordsr  exist  (fables  5  and  6)  red  slice  analytic 
evaluation  of  effects  of  system  uncertainties . 

tl**,**,™ 

Tor  aecond-order  system,  5  i  0.65  for 
n  ■  2.  (Ref.  57)  Selectivity  bcccmes 
greater  u  n  Increases,  and  criterion  alto 
become*  wore  complex  to  aaaeal. 

Analytical  expression  be  case*  more  difficult  •*  u  Increases,  vlr: 

X”  tvWd. .  (-»» 

Standard  forec  and  associated  transient  responses  have  been 
evolved  for  n  •  1,  2,  5  using  cceputers.  (Ref.  70) 

General  (1  «»,  nonmintdied)  fora*  presented  In  fable  5. 

Standard  fores  In  fable  6.  Generally  not  worth  the  extra 

computing  effort  over^°° ta?(t)dt. 

mm 

Held*  (  •  0.667  for  secood-orler  ays  tea 
fcr  c—  co.  (Ref.  75)  Haa  not  yet  been 
evaluated  for  higher  oroar  system. 

Cxtrenn  vrlue  can  be  found  analytically  (fief.  75),  *1  the  ugh 
even  sixple  fores  are  fairly  coagdex.  Selection  of  a  la 
arei»rary. 

It  another  parameter,  like 

(torque)^dt  »  c<xwt. 

Ralls  the  ready  appllc ibUity  teat  for  the 
constant  coefficient  ayrtese  considered 
here.  Choice  of  sl,  1*.  arbitrary,  ashing  the 
general  fora  of  ihlo  neusure  difficult  to 
assess  In  the  present  teres  of  reference- 

A  typical  crlt-rim  fore  suitable  for  solution  via  dynmic 
programing  technique*.  This  one  has  been  selected  frost  a  large 
number  of  relatively  unevaluated  "ganeralltatlona"  proposed 
recently  (iWf.  7,  8,  9»  51 ,  *6,  *7)  because  it  reduces  In 

Halting  cars*  to  criteria  that  have  writ  for  staple  constant 
coefficient  system, 

«.g>,  alnlaai  r c2<5t  vhlle  ( torque }2dt  -  const  This 

select*  (  •  0.7  for  a  second-order  aystc*.  (Ref.  63)  Qae  not 
been  evaluated  for  higher-order  system.  Will  be  cosh  more 
fsailtar  and  important  as  dynamic  protreealt*  application* 
ecjaM. 

2 

Frobaoly  fails  ready  applicability  teat. 

Choice  of  t  la  arbitrary,  asking  the  general 
fcn>  of  this  ssuirt  lil-defined  and 
difficult  to  assess. 

r[s,  t,  lj]  la  a  general  functional  if  •rror,  tlsn,  and  system 
parameters,  Ag.  p(t)  is  the  probability  that  the  output  will  be 
used.  Proposed  a*  an  *11 -encompassing  criterion. 

Iks  value 

ASSOCIATED  nSRCKS*  .*CAS0SR3 

Sot  ywt  evaluated.  A  logical  fora  for 
ays  tens  with  finite  position  error.  Choice 
of  T0  la  arbitrary. 

A  clever  use  of  Uguerre  function#  In  evaluating  the  integral 
appear*  In  Ref  3,  which  also  cox&lders  the  cece  with  an 
lapulae  Input.  Kay  have  application  to  special  control  problsm. 

m-nta.*-® 

lot  yet  evaluated,  lock  of  fined  value  for 
?„  sake*  this  aeaaure  vague  and  lodeflnlta. 

Proposed  «  a  raeerallted  oerfominc*  aeaaure  Tor  transient 

inputs.  Reduces  when  70  •  0  to  e2(t)dt. 

X”  '.<«>« 5 

X° 

reSiit  ptedstendned 
lint 

Unsatisfactory  response  for  higher  order 
system.  (Ref.  » 

7 imply  re le ted  to  irput-outjwt  transfer  function  In 

asm  way  ax  error  coefficient#  are  related  to  tlm-welghted 
lapel?*  response  of  error  (see  Chapter  V). 

2 


is 


1 


(29) 


£lS 1  .  _ 

a<5)  i,2  *258+1 

For  systems  with  a  nor.ze:-o  cteady-ctctc  error  following  a  step  Input,  indicial 
error  measure  aic  uol  explicable,  as  novtd  aoove.  Hcve/tr,  it  Is  easy  to  derive 
associated  '“•“C'nied  response  measures  yielding  finite  criteria  values;  four 
such  Erasures  are  appended  to  Table  IV. 

The  aeasures  will  now  be  examined  with  regard  their  validity,  selectivity, 
and  ease  of  application. 

An  noted  in  Table  TV,  tk,  jte,  and  lT^  are  all  invalidated  by  thejr 
inability  to  discriminate  between  responses  which  are  grod  in  the  sense  that  the 
absolute  error  is  snail  or  d-cay3  rapidly,  and  oscillatory  responses  which  are 
lightly  damped,  and  in  which  the  large  positive  ana  -.egatixe  errors  approximately 
cancel-  Similar  drawbacks  apply  to  the  associated  ixpulaive  response  aeasures 
/^(tjdt,  tnc(t)dt  Therefore,  tn*ri»e  aeasures  will  not  be  discussed  fur¬ 
ther  in  tnic  chapter,  except  insofar  as  they  represent  limiting  value:*  of  otner 
measures  (e.g.,  ITAE  -  IE  for  nenoversbooting  responses),  ifoit  of  the  rcaair.-der 
of  this  chapter  will  be  devoted  to  a  stuay  of  the  Was  ares  IE2,  ITE2,  IT^E  ,  IAE, 
ITAE,  and  IT^AE.  Particular  attention  Yiil  be  gl  -en  to  the  analytic  evaluation 
of  *hese  measures,  net  only  bv'o.ase  of  the  insight  tnat  analytic  expressions 
juoviue,  but  also  because  they  enable  a  cnect  ta  be  male  on  published  values 
obtained  by  mechaniiat’.on  of  analog  computer  re-poncts.  As  will  be  shown, 
several  error's  have  been  detected. 

DERIVATION  OF  IE2,,  ITE2,  AND  IT^E2 

Tables  giving  IE2  for  nonunit  numerator  systems  of  first-  through  as'.with- 
order  in  terms  of  numerator  and  denominator  -*oefficie"**  are  given  in  Ref.  V., 
Appendix  E  of  Ref.  58  extends  the»*e  tables  through  tenth-order,  and  corrects  an 
error  in  the  scventh-o;der  integral  vf  Re«*  44.  These  literal  expressions  are 
lengthy,  ftiiu  n.-f .  bb  notes  that  tb.  integrals  con  be  expressed  *wre  uuapactly 
in  terns  of  Huruitz  determinants,  a  point  which  is  discussed  further  in  Ref,  6 
and  4p.  Tables  for  ITE2  (*nd  IS2'  are  $Jve:*  ly  Jfe-ieotl  (Ref.  79)  for  systems 
ol  first-  through  fourth-order-  The  topic  also  been  studied  by  Knothe  (Ref. 
48)  and  Stone  (Ref.  70),  using  a  ecwvhat  different  approach  tlwn  t^e  previous 
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references.  Stone's  procedure,  and  that  of  Veatcott,  is  outlined  very  briefly 
below. 


She  procedure  is  aiaple,  in  principle,  once  the  literal  fur*,  for  Sq  y\  has  been 
obtalnM.  However,  in  practice,  for  systems  abftv*  *v,>  ♦*>«<>  iu»^i 

fcm  becomes  exceedingly  long.  f Stone  presents  such  a  fora  for  a  nonunit  nuser- 
otor  fourtb-ordv-i  system  which  uccupiec  nine  pages.)  Ccac  3 ir^lifi cation  any  be 
introduced  into  the  differentiation  required  to  obtain  IT^E^  by  the  a  priori 
neglect  of  terns  ir*  and  above,  which  Vardsh  when  s  is  allowed  to  tend  to  zero 
in  the  second  derivative.  iWertr-eiess,  Vestcott1*  procedure  seems  to  be  cc«at- 
vhat  briefer,  end  is  preferred  for  the  purpose  of  detaining  IT^and 


* 


.•H 


A 


»*  f. 


I  * 


V9SW5  *  i'oli.'Wing  Pef.  ^  bud  58,  eo-vioys  Parse va 1 ’s  theo-ei.  to  express 
as 

IE2  -  J”  [e(t)]£  dl  .  „^0  gJJ  j^+Ja>  EC*)E<°  -  =>)d«  <»> 


which,  for  u  hlnble  »jste=.  beccce* 


9  1  f*-100 

XI?  •  ggj  E(B)«(-8)d* 


When  B(s)  1#,  in  sdditirn,  ft  rutlo  of  rational  polynoniaVs,  Eq  36  la  •  syaaetrl- 
cal  rational  function  of  the  poles  of  E(s) ,  being  the  sua  of  the  residues  at 
Usee  pole*.  *u*i  can  therefore  be  erp rea^A  In  tern-  of  the  coefficients  of 
£(s)  only,  The  procedure  is  veil- illustrated  by  considering  the  simplest  case: 


which  has  ft  single  pole  at  s  »  enclosed  by  the  contour  integration  which 

is  an  arbitrarily  large  *c?ai circle  in  the  left-half  plane.  The  residue  at  this 
polo  is  equal  to  TE2,  and  is  siren  by 

ijS  .  “a  [*0*  *  *'  .  1  .  j£_  (w) 

,-jSlL  *0  (Bo*  ♦  -  «1>  J  2v>, 

The  procedure  for  obtaining  ITS^  U  indicated  b7  Eq  ”• 

*  f  <•[«(*■) J"  at  - -c^osl i7$fc  m  -  «)d*  i*o> 


It  ic  not  possible  to  take  the  o,  _►  o  Uniting  process  under  the  integral  sign 
without  first  perform- ng  toa  differentiation.  To  obtain  a  ryswetrieal  function, 


-c/  • 


assumed  to  be  real  and  jCPitire,  the  contour  is  teken  as  the  mid  line  of  the 
ulsplacx-went  of  the  E(s)  and  E(-s)  functions,  completed  by  n  large  semicircle 
occupying  the  left  £slf(  and  pa Tk.  of  tb*  right  half,  ccrylex  pisne.  Ibis  ie 
equivalent  to  sinply  changing  the  variable  a  to  s  +  (oj/2),  aod  treating  the 
contour  Integral  as  an  Ordinary  Integral.  Equation  J6  theu  beccaea 

*~+  ,»0> 


-  &  Si  ■(•♦tKV  ■5)d* 


By  substituting  o  »  Oj/2, 

=*  * -**[*C>*~-H 

This  nay  now  be  solved  by  naklng  use  of  the  e£  tables.  For  eraqvle,  taking 
the  wispiest  Cate, 

E!‘>  -  w 

the  intrgrana  becomes 

_ 4 _ 

(»0*  *  >0'  +  «|)(-«o»  »  «oO  ♦  a,) 

The  group  a1  ♦  e$o  now  correspond*  to  a^  In  Eq  38;  Caking  the  appropriate  substi¬ 
tution  in  Eq  VO  yields 

_ 2  I  llJ  S  /  4  \  do 

‘  ■  2  53  [XqI.w  -  .,}/  ■  £2  (1,,) 

By  means  of  Uila  procedure,  Weatcott  obtains  literal  forms  for  ICS'*-  The  process 
ha*  been  extended  by  a  further  -d/&?  operation  to  obtain  IT^E?  in  the  present 


Table  V  presents  the  resulting  literal  fonts  for  TE?,  IIP?,  and  Tjri?  for 
systems  of  first-  through  fourth-order.  *h«  iI5?*  results  hsve  been  ralerived 
independently  of  VestcOtfs  volt,  srd  check  his  values  exactly.  Jto  genera) 


urn  v 

liter',  l  raoe  cr  nn&auis  respired  rot  ie2,  ite2  led  e^e2 
TOR  SIS?S>S  <y  FUST-  TEROWH  TOURTH-CSDES 


rr/  v,f 

'  i.  r  -sj 


*>  *  4* 


*s4  •  v*?  -  =v.i  ♦  ^  4 


4<*S*3  -  *■•»>  *  *3*J<4?  -  *V2>  •  *0H<4  -  «|’>>  ♦  ^  *5  <*>*»  -  vp 


4  l  '<*i*a*5  *  *o*5  -  »f«J 

^  («1»;  •  1j*p  "  S~  '*?  *  *C*»>| 

SC*rt*J  •  «0*J  -  4*J  J 

(4  •  *i«j  -  *  «jm*?  -  «j«j)  ♦  •?(*!  ■  *mj>  •  5  -  »*p] 

S<*|*3*J  -  «e»j  *  •?[»)* 


TABLt  V  (Continued) 
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Xa 


-  5^ 


k{< 


«C«1<  -  -  Sfibh 


■t«f < s>4«S  ♦  ^44*?! 


>»« 

41  *  2^)£*o^*5  *  *+3  *  *|4> 

2j  •  S/M*)  •  •  *t»*> 

*5 

*5  *  •  W*1?  -  >Vj>  •  V,^ 

«j  .  •  «»i*J>  ♦  Vc*-?  -  ®W  •  •  5-o"i  •  *jW 

>j  -  ^  .  Wp  *  *«g*tO*<ftt  *•>«»•  *!«0> 


«i  •  *"o<»i*e*5  -  vS  -  £?•»> 
k  *  *V>-’’a  ♦  *1*3  -  *M*> 

o  •  *  -*i  •  -  *4»> 

k  •  •  «u*3<*f  -  =Vi’  '  *o*(<*i  -  “>«;>  *  5  •  *o*J> 

k  •  »€(•! - •,»> - »mh; •  •§[‘4 - «i*> 

•  ■<*[»‘i*t  •  Hs*j  •  9  2  Vj  *  ‘i>»  5  ■  *  »W*?  -  *wV  * 

5.  -  xjtfci*  ♦  -  3VP  -  s  (4  ‘ 

•  w|,c*f  •  '**  • * !f  <4  •  “i*P  • 

•  <Wi 
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literal  values  of  IT^E^  have  pievicusly  teen  published;  hence.  It  is  not  post«ble 
to  have  a  -voplctc  independent  check  on  this  measure,  olttaish,  as  will  be  noted 
ir  tr.c  following  section,  n  partial  check  does  exist. 

KVALUAJ10K  C?  1*?,  m?,  AND  IT2*?  AS  FKHFCWWKTB  MSASTKkS 

Indicia!  responses  of  normalized  unit- numerator  mini mas  IE2,  IT!?,  IT^x?, 
and  IT^E2  systems  are  given  Is  Ref.  70.  These  arc  reproduced  is  <  lg-  ?5  for 
ease  of  reference,  ana  corresponaing  closed-loop  Bode  diagrams  are  siren  in 
Fig-  !*».  The  transfer  function  coefficients  that  minimis:  these  normalized  IE2, 
IT!?,  etc.,  performance  measure*  are  called  standard  fores,  i.e.,  the  partial 
derivative  of  th<*  normalized  performance  measure  with  respect  to  each  coefficient 
of  th*  eloivl-lccp  transfer  function  ia  zero  vhen  the  transfer  function  nas  the 
appropriate  stancard  fore.  Table  vi  present t.  such  standard  form,  for  a  variety 
«f  per f^rrance  criteria,  including  minimum  IE2,  IT!?,  and  IT^!?-  Those  have 
been  ootalnec  by  dlffereu.iating  tbs  general  analytic  .-inressiens  cf  Table  V 
with  respect  to  the  transfer  function  coefficients.  The  results  "crec  with  these 
of  Ref.  70  tihlch  were  obtained  by  a  limit*!  computer  programed  for  iterative 
plnli^LsaUon.  (This  provide*!  a  partial  cheer  os  tr.c  accuracy  of  the  It2!?  gen¬ 
eral  forms.)  Reference  70  eJao  presents  standard  furo.  and  xndici&l  responses 
for  optimal  VthP  syuteaa.  These  wlyc  not  been  cnecaed,  tut  are  inclulea  in 
Table  VI  to  ensure  completeness.  The  anal^lc  expressions  for  IT^e2  for  sysUm* 
of  higher  tha.5  ccccnd-order  are  very  cotsplicaW,  and  IT^?  does  sot  appear  to 
offer  advantages  over  P'S?  or  IT2!?  that  would  he  sufficient  compensation  for  the 
computing  efforts  involved  in  practical  optimization  calculations. 

Tible  VI  also  lists  the  open-loop  transfer  function  coefficient*  associated 
with  the  standard  fores.  As  with  the  closed-loop  coefficients,  the  sensitivity 
of  t1*  normalized  performance  measure  to  5 sale,  changes  in  theoe  coefficients 
e*x>ut  the  specified  values  is  se*«.  Hw.  uc^erer,  teat  the  unno realized  per- 
f crrAT.ce  measure  is  effected  by  szall  changes  in  open-loop  coefficient*  from  the 
standard  forms. 

Trie  ciosea-iecp  ( Juj)  co,ic  aiagrams  ror  tne  xz?,  !£?,  n?s?,  ana  iT*z? 
standard  form  arc  given  in  Fig.  I5*-  (5  templates  correct  to  one  significant 

figure  vero  used  to  co*istruct  these  diagrams.)  By  operations  cn  a  Jlichols  chart, 
or  by  direct  factorisation  of  the  lower-order  systems,  the  associated  open-loop 
•:a '■ceteris tics  of  phase  margin,  crossover  fren.je r.cy,  etc.,  have  been  found. 


«0i 


I 


These  "re  presented  ii.  Table  VII,  together  wIth  th*  corresponding  Y»hAs  for 
ninlwes  I7AE  systems  and  other  standard  fores.  Table  VII  and  Fig*  1.1'  indicate 
that  the  standard  forms  having  the  smoothest  indicial  response#  posses  phase 
margins  between  60  and  TO  degrees. 

For  a  no  real!  zed  second- order  unit-numerator  xj*lv.  1  wiving  a  transfer  func¬ 
tion  aescriW  by  Eq  29,  «  -  7'^%  Toe  general  expressions  for  IE^, 

ITS2,  and  n^E2  cirpllfy  to 


IE2 

„  1 
-  ?  *  q 

P*> 

. - O 

lit*” 

.  i’.-4 

(*5) 

6?‘ 

r^E2 

- 

e?5  B? 

Pti.lt 
-  ^  - 

fi*6) 

These  measures  have  been  evaluated 

and  graphed  in 

n«.  i£. 

Tne  grapns  of  Ref.  37,  which  were  obtained  by  mechanization  of  analog  ;ocs- 
puter  responses,  are  also  shown  in  Fig.  l6  for  comparative  purposes.  The  dis¬ 
crepancies  in  these  •  .irves  are  likely  to  be  attributable  to 

1 .  A  scaling  error  of  a  factor  of  1 0  for  ITE2 

2.  A  scaling  error  of  approxirately  l6  for  IT^E^ 

3*  Amplifier  drift  fer  £  >  0-9  for  all  three  measures 
The  conclusions  of  Ref.  37  regarding  the  value  of  these  measures  as  criteria  for 
normalized  second-order  systems  i**ed  not  be  changed  by  the  correction  of  ibese 
errors.  IS2  is  minimized  by  £  •  0-5,  which  is  not  obviously  undesirable,  although 
rather  higher  damping  ratios  (5  «  0-7)  are  usually  preferred  for  step  responses. 

The  strongest  objection  to  IE2  ss  a  performance  measure  for  second' order 
systems  is  its  lack  of  selectivity.  Varying  *  free  0.2  to  1 .3  raises  IE2  frea 
its  minimus*  value  of  1  to  only  *.*».  By  contrast,  ITE2,  vaich  has  a  niniaa  value 
of  0.9  at  £  ■  0.67,  increases  50  percent  (to  1.35)  at  £  »  0.51  *md  5  -  0.5 .  Per 
the  normalized  third-order  system  investigated  in  Ref.  57  with  the  transfer  func- 
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IE2  tcccsscs  efcn  Ires  selective.  As  ch or  ir.  Fig.  11  of  Ref.  37,  for  b  ranging 
iroa  1  .25  to  3.0,  c  can  be  vrrled  froa  1  -5  to  3«0  «itii  less  than  10  percent  re¬ 
sulting  var:atlc:.  In  I£^.  The  Indicia!  responses  of  Fig.  15  Indicate  that  oint- 
sisotjor.  of  II*-  selects  oscillatory  responses  lor  highcr-onier  systems. 

An  advantage  of  IE^  is  tnat.  its  c..tl/lic  expressions  «r*.  simpler  than  those 
Tor  fTF2  and  TT^E2.  The  last  seas  ire  is  highly  f  elective,  but  it  is  doubt^il 
vhethe;  the.  added  ccnpiicot<ou  of  the  general  fores  is  wortu  the  extra  coeputa- 
tio  >al  effort  required,  csrpared  to  I7E2. 

Having  stated  the  advantages  nr.d  disadvantages  cf  12- ,  ITL2,  and  IT^'E2. 
sect  alternative  criteria  can  nov  be  assessed  to  arrive  at  a  unique  selection, 
if  such  a  selection  can  be  rjLdc.  7t>>  next  section  of  this  cheater  is  concerned 
with  TAR,  I7AS,  and  hut  it  '«uli  be  liiolcading  to  close  tue  discussion  of 

IE2,  ITS2,  etc.,  h  '^iUT  viigrecalori  relating  to  the  suitability  of  these 

criteria  for  statistical  inputs. 

A  strong  argument  in  favor  of  IE2  is  t>.at,  for  statistical  inputs,  the  cor¬ 
responding  measure  e^(t)  •  ^  [e(t)]  ju  hxs  been  thoroughly  investi¬ 

gated,  and  foms  the  foundation  '‘or  an  extensive  literature,  Including  *+ny  cf 
tnc  cost  valuable  and  widely  used  contributions  to  optimisation  theory. 

Statistical  inputs  fail  outslae  the  scope  of  present  report.  However, 
it  vculd  oc  unrealistic  to  ignore  the  fact  that  r light  control  systees  ere  iuo- 
Jeeted  to  both  statistical  and  dateraiatst i *  inputs,  and  it  would  ee  unfortusetc 
if  the  already  considerable  gap  between  mm  dee*,  and  deterministic  setooas  of 
analysis  vat  widened  ..y  dcsrrda  for  different  critcila  for  each  type  cf  Input. 

It  is  suggested  tbrt  the  possibility  of  developing  a  <*-x  eatable  criterion  Right 
well  fern  the  subject  of  a  future  investigation.  However,  for  the  remainder  of 
this  repoc*,  statistical  considerations  w ill  w  disregarded,  nnd  assessments  of 
Indicia!  error  oeasorcs  will  be  nsec  with  respect  only  to  their  suitability  for 
the  tosh  of  measuring  system  performance  fol'ovirg  a  3ti*p  input 

rEHivATicn  cr  nb,  itaS,  ah^  h^ai: 

The  1AR,  flAE,  and  IT^AS  measures  have  beer,  investigated  in  ?cf .  3?,  36,  and 
X>.  Reference  37  includes  an  extensive  study  of  the  properties  cf  ITAE  obtained 
by  rscchanlratioe  tf  a-aloc  computer  responses  for  iwit-nirrcrator  sero-position- 
error  syatesA  of  first-  thmueh  eighth-order,  first-order  nunc rater  tvre -velocity- 


.  -i  n  v '  .  -ow.«-  ,  u-.wgn  sixtn  -^rQer,  ood  zcro-acceleration-error  systems 

of  third-  through  sixth-order.  indi^ial  rtspoi-scs  and  standard  fores  for  sach  of 
those  oysters  an  ^resented.  Ikr*.-  fjtta  are  also  given  on  IAE  and  IT^AS,  but  TVAE 
Is  enponslzcd  because  it  If  *»**icct've  chan  i\S,  and  is  easier  to  apply  than 
IT^AE-  (ITAE  cat.  L*  easily  cechaj'zca  ana  scaled  on  an  analog  computer,  c.g., 
using  the  simple  circuit  described  Jo  Ref.  Vf.) 

The  present  report  continues  th*  l ivestigetion  of  these  measures.  Exact 
'Viaiytic  expressions  for  IAS,  ITAE,  ar.l  IT^AE  or  *econl-c-r»?cr  »ystea s  ere  pre¬ 
sented,  the  res'tltc  are  employed  to  check  those  of  Ref.  37.  As  vill  be  sfcovn, 
•ore  ^rrs'.r*  arc  d* tested.  An  analy'ic  method  rcr  obtaining  IAE,  TTAE,  and  ET'AE 
for  hicn^r-^'dcr  systems  is  also  -h«Mr.f^  in  addition,  Chapter  *  j  dvocxiceo  - 
procedure  for  the  approximate  calculetu  .*s  cf  ITAE  for  near -opt  lain  systems. 

The  procedure  for  calciiatirg  IAF,  ITA2,  and  IT^AE  for  a  see<hri-o«»~r  y^- 
tes.  Js  detailed  in  Appendix  ’ ,  and  brfafTy  s^usiarised  bcloo. 

The  starting  point  is  tre  calculation  of  the  Laplace  transfom  cf  the  ab¬ 
solute  error  tioc  Mctc**y,  A  (?;  Once  in':  bar  been  obtained,  IAE,  ITAE,  and 
IT^AE  ere  r-edtty  cotaincd  by  peans  0/  ttv*  final  value  theorem,  as  indicated 
bolov. 


IAE  * 

to 

t|-®J0 

i«(t)!dt  ^ 

to  eA<c> 

»->c  a 

(W)) 

ITAE  - 

iu  /*< 
•iW. 

t|e(t.:dt 

to  1  •dElt*I 

*  S  -►  0  s  ic 

(*9) 

IT^AS  > 

to  /-1' 

to  1  J'tA(c) 

*— *  »  4,2 

(50) 

The  method  for  obtaining  £a(s)  anMytiOfiliy  is  an  extension  cf  the  tech¬ 
nique  u*ad  '“or  obtaining  the  I*p?act  vr^sforn  of  »  rectified  sine  wave-  Ae 
error  time  history  lie  damped  sine  wave  v'.tu  a  phase  lag 
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From  the  ..ransforn  tables  of  p.-f.  lo, 

4.!^!r.^l  .  *  (-0'n  <;j> 

As  shsvr.  in  App*>odix  A,  tic  damping  factor  la  Eq  >»  can  be  accounted  for 

by  replacing  s  with  (  ♦  £0^,  to  obtain  <at|.  The  phase  lag  intro¬ 

duces  27Zt  ecri  lira* ions.  bue  the  final  result  Is  reasonably  conoact. 

♦^^/co  p  a(j  4.  ^Ow)  *]  s  ♦  2^ct, 

£a(»)  -  - -5 - 3  COth - *  -  V-  :  {» 

(s  ♦  ♦  o*  «■  J  (s  ♦  ta.)  ♦  co* 


-  ct  V »  -  S2 

*Ote  that  the  la**  tens  in  5^  is  simply  E(«). 

tts  procedure  In  obtaining  an  analytical  expression  Tor  |e(t)|,  or  Sav*)# 
for  yefic*al  higher-order  systeas  is  jore  complicated.  For  these  systems,  the 
the  7Civ  trotwlMS  of  .he  error  tins  history  are  not  equally  space!,  and  Eq  55 
d'-es  cn>  apply.  '*t/-  procedure  used  in  Appendix  B  for  a  third-order  system  having 
o  c  cep  leu  pair  c.‘  roots  *rvs  to  rspreuent  the  error  response  by  a  Fourler-1  th¬ 
eories,  the  coefficients  Of  which  are  t lac-dependent  The  final  result  1*  quite 
c  cep  Heated.  It  pruvidrs  a  basis  for  further  theoretical  studies,  and  for  check¬ 
ing  ef>cria^rtally  obtained  IAS,  etc.,  but  seems  too  involved  for  routine  opti¬ 
mization  calculations. 


por  a  normalized  unit- numerator  secowi-oroer  syz; 
IAE  and  17  AS  are  given  by 


•cm  described  by  Eq  2y,  toe 


hfe}"]-1 


end 


•rn?£i  w*w»*»ssi<ma,  tq  a**«o  or  Api-e*vlix  A  Jwhicb  It  s^.*..  wet,  TlfAEi  «a«s 
scnted  in  Fig.  17  a«d  I 8,  and  art?  compared  with  inc  vale**  obtained  ty  Grahaa 
and  Iiitnrop  i*i  Ref.  $7-  It  Ip  evident  that  *  scaling  error  of  tvo  sxists  in  the 
•«nh«n  rrl  lytlliro*.-  yntyii  for  Tt^KZ-  Agrees nt  on  IA£  is  excel lone;  sock  rslncr 
discrepancies  exist  in  ITAE.  Klsl&lzatlou  occurs  at  £  »  0-76  instead  of  £  *  0.7, 
but  this  error  is  insignificant* 

The  conclusions  of  Crahaa  ar,i  Iathrop  regarding  the  relative  nerlts  of  these 
criteria  arc  not  altered  by  the  present  study.  As  vould  be  anticipated,  IAE, 

ITAE,  and  r^AE  exhibit  sicilxr  characteristics  to  IE2,  ITS2,  and  TT2^,  respec¬ 
tively.  IAE  is  uwselectivc,  and  IT^AE  is  tco  cccplicutedi  therefore,  TTAE  is 
preferred-  ITAE  has  been  explored  t«orc  thoroughly  than  any  other  indicial  error 
seasure,  snd  the  considerable  background  of  knowledge  established  by  Ref.  >7  md 
?8  pemits  Ihir  ncasure  to  be  used  *'ith  confidence.  For  each  case  investigated 
in  Ref.  yj,  rdnlnun  ITAE  yielded  "good"  responses,  and  tne  criterion  regained 
highly  selective  for  cystess  of  third-,  fourth-,  and  fifth-order. 

Co-parcd  with  HE2,  ITAE  has  the  advantages  of  giving  lest  oscillatory  in¬ 
dicial  standard  foms  and  of  bei  g  sore  thoroughly  explored,  but  it  is  the  acre 
difficult  of  the  two  censures  to  express  analytically.  The  choice  between  these 
criteria  cult  depend  on  the  circumstances  of  the  individual  cases  in  which  they 
arc  to  be  applied.  Possibly  ITE2  would  be  aore  convenient  for  an  analytical  in¬ 
vestigation,  and  TTAE  best  for  eptinizatien  using  analog  computers,  because  of 
the  ease  with  which  is  can  be  ccchanized  and  scaled. 

Closed-loop  Sole  diagrams  of  the  ITAE  standard  fora*  listed  in  Table  VI  are 
given  in  Fig.  1*».  Associated  indicial  responses  are  given  in  Fig.  15*  Table  Vin 
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Comparison  cf  Analytical  and  Expert,  non  tally  Determined 
1AE  and  1TA1  fop  Unit-Nuseratar  8ec3f  <l-Order  System 


TABLE  YUI 

ECJICIAL  EESPCHBS  CKARK7ERIOTICS  OP  CPUKU  StOTBC 


' 

Delay  Tiae 

Else  Tiae 

Time  to 
Peak 

Peak 

Orerahoot 

XquirjLlent 

Tiae 

Constant 

Settling 

Time 

Second-Order 

|| 

Minimum  Settling 
Time 

l.*3 

2.2* 

*.*0 

0.066 

n 

2.96 

1-31 

8.5*> 

5.03 

0.023 

mm 

3.67 

J? «C6>  *6 

1.38 

1.65 

3-68 

0.165 

el 

3.36 

eft)  2dt 

1.6 

1.86 

3-93 

0.096 

1.67 

5-86 

^“t2  e(t)  2i)t 

1 M 

2.00 

6.85 

0.068 

1.78 

5-03 

jj“t3  eft)  2it 

1.5 

2.21 

4.6 

0.061 

l.8i 

3.05 

Third -Order 

Minimus  Settling 
Time 

2.02 

8.03 

*•32 

0.031 

2.32 

3.55 

£°°t|e(t)|4t 

2.1! 

83 

6.67 

0.029 

2.** 

3-38 

lW  eft) 

1.8J 

1  -76 

3.62 

0.071 

2.10 

C.7I 

JTrot  eft)  2dt 

!■!* 

1.95 

>•90 

\06i 

2.5»i 

6.98 

X^t2  «M  ®4t 

2.0 

2.32 

6.16 

0.065 

2.30 

6.82 

«<*> 

2.0J 

2.20 

6.56 

0.038 

8.37 

3.61 

Minimis  Settling 
Time 

2.56 

8.17 

6.78 

0.0? 

8.88 

3-8 

jT^tleftJUt 

8.62 

8.Vf 

5-66 

0.027 

3-C8 

6.25 

j£“  «l»>  *4* 

a.6o 

1.8J 

7.75 

0.161 

8.68 

10.32 

jf’t  eft) 

2.1>9 

1-5* 

8.02 

0.06C 

8.76 

8.91 

£"t2  eft)  «4t 

2.58 

2.16 

|  | 

0.056 

8.86 

5-0 

£®tJ  eft)  3it 

2.62 

-  . 

2.23 

o 

0.060 

_!f_ 

3-93 

60 


suETArizet  the  indicial  response  characteristics  of  all  the  standard  foiw  of 
Table  VI.  Standard  fora  root  locations  are  shovn  in  Fig.  19,  20,  «M  21. 

onu?.  IKDIC1AL  f-Kx.C*'. 

The  measures  discussed  in  the  previous  section*  of  this  chapter  have  all 
been  Integrated  functions  of  error.  It  1*  possible  to  generalise  these  Manures 
by  Including  functions  of  the  tlse  derivatives  of  error.  Thus,  several  recent 
references  (e.g.,  Ref.  79  anl  d0)  have  proposed  optimization  procedures  based  on 
the  fell  owing  very  general  measure: 

f  fte(t).  t»  7t]?it)4t 

*  U 

Here  F[c(t),  t,  J  is  a  general  function  of  error,  tlse,  and  system  parameters 
7 it  and  pi  is  the  probability  that  the  output  *111  be  used.  The  generality  of 
this  measure  Is  both  its  advantage  and  its  drawback.  It  cannot  be  put  into  con¬ 
crete,  usable  for*  without  making  a  mafcer  of  arbitrary  choices  to  arrive  at  F. 

Consideration  of  the  remaining  measures  listed  in  Table  IV  illustrates  this 
point.  The  criterion 


yields,  for  a  norrxlited  second-order  unit-nuaerator  ays  ten,  £  •  0.866  for 
a  -»>  a>  and  0-5  for  c  0.  (As  shown  in  Ref.  75,  the  £  selected  by  the  crite¬ 
rion  can  fall  out-  .de  the  »-ange  0.5  to  0.866  for  a  negative  c.)  Ho  guarantee  is 
available  that  *n  a  suitable  for  a  second-order  systea  will  yield  acceptable  re¬ 
sponses  when  the  criterion  is  applied  to  a  systea  of  higher  order.  However, 
the  criterion  cannot  be  reflected  on  this  basis  alone  without  further  investigation, 
particularly  as  its  randoa  analog  should  be  easier  to  use  with  statistical  inputs 
than  the  randoa  analog  of  tine-weighted  criteria. 

Arbitrary  constant*  also  appear  in  the  quadratic  fore  associated  response 
measures  listed  at  the  bottcc  of  Table  IV,  and  ir.  the  very  general  neasure 


"Igure  19.  Pole  locations  of  Unit- Kuxera tor  Second-Order  Standard  ?or*s 


Figure  21 .  Pole  Locations  of  Unit-Hunerato/  Fourtn-O.iier  Forms 
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Here  the  choice  of  aj,  ag,  . ..,  an  is  arbitrary,  so  that  any  desired  weighting 
of  the  elements  of  the  criterion  may  be  used.  Unfortunately,  little  guidance  ia 
available  on  the  choice  of  a, ,  ag,  •  ,  an  and  considerable  effort  aay  be  required 

if  it  is  necessary  to  evaluate  the  integrals.  (However,  this  measure  has  usually 
been  proposed  in  conjunction  with  dynaaic  programing  procedures  employing  digital 
computers.  Kor  such  applications,  the  latter  disadvantage  would  be  insignif leant.) 

Hew  measures  are  constantly  being  added  to  the  already  vast  literature  on 
this  subject..*  It  Is  perhaps  worthwhile  adding  a  plea  for  caution  in  accepting 
these  measures  as  valid  criteria.  Tt  ia  all  too  *asy  to  present  a  "rev"  per¬ 
formance  criterion  which  is  valid,  selective,  and  even  readily  applicable  fer 
scoe  particular  systen.  The  task  of  formulating  a  criterion  which  retains  these 
merits  ever  the  range  of  systems  encountered  in  flight  control  design  is  much 
more  difficult.  Of  the  indicial  error  measures  dlaeusscd  in  this  chapter,  ITA2 
stands  out  because  of  the  thoroughness  with  which  Jt  has  be«n  investigated.  The 
large  number  of  examples  given  in  the  studies  or  Graham  and  Lathrop  lend  credence 
to  their  assertion  that  the  ITAE  criterion  has  exceptional  meric. 


*Thl&  point  is  well-illustrated  b,  I5f  f.  42,  which  was  received  while  the  present 
report  was  being  typed.  In  addition  to  IAE,  IE?,  TTAE,  ITE,  ITE?,  and  li?AE, 

Ref.  42  investigates  the  following  measures  for  normalized  unit-numerator  second- 
end  third-order  systems: 


f  VUWl  dt,  f  Vt|*(t)l  dt,  f  t»/|e(t)|  dt, 

•'o  ■'o 


r  ^vrwt  dt.  r . 

•*o  j 0  1.  [e(t)]* 


Both  digital  and  analog  computers  were  employed  to  generate  the  measures.  The 
results  are  generally  in  agreement  with  those  of  Ref.  70  and  this  report,  except 
for  errors  (in  Ref.  42)  of  approximately  10  percent  on  acne  standard  form 
coefficients. 
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OTHER  STAIIChHD  FOhw; 


lac  standard  s»»r»«  discussed  *_*?rli«»r  In  this  chapter  have  all  been  cptioal 
L*“  Uk  sense  that  they  '•"•Ut  in  Binlfcue  voices  of  particular  judicial  error 
r  •i*ir«ria  could  Lc  used  for  op lini nation,  ana  riex.  jd  pi'.^c.its 
tat  standard  forrrs  associated  with  ninieua  settling  tine  for  systens  of  second- 
thruurl  t  ighth-ord"r.  The  cccona-,  tnird-,  aj»d  fourth -oraer  tons  ana  associated 
transient  responses  art  r^ro^nwi  in  7**VI»*  VT  *na  fip.  15.  Figure  1 4  illustrates 
the  lode  (  Vo)  diagrams  re:  these  systems.  Table  VI  also  lists  binualal  and 
Sullcrvorth  stPnutud  itms.  The  oincclal  fores  are  not  optiii&l  in  the  some  sense 
ac  toe  ctar-dar-l  forts  considered  previously,  being  generated  by  the  following  cx- 
pn*“.  .Son: 


Sfs} 


(0  *  -V" 


(sv) 


»•  is  the  aystec  order. 

The  Buttervortu  foms  arc  cliaractcrizcd  bj  the  lact  that  their  nolcc  arc 
equally  spaced  around  a  circle  of  radius  ft©  in  tut  cccplcx  plane.  The 
frcqutjjty-dccain  representation  of  th*  resulting  t-rarafrr  function  iv,  illustrated 
in  Fig.  1’*  for  second-,  tim’d-,  and  fo»i»th-ordcr  syotece,  the  Indlclal  responses 
arc  given  in  Fig.  Vj.  Buttervorth  fores  arc  maximally  flat,  i.c.,  at  Jee  0,  the 
first  n  derivative.,  ui  cn»-  Mqplitude-freqnor.cy  diagram  ere  zero,  but  their  tran¬ 
sient  responses  have  increasingly  large  overshoots  «s  the  order  of  the  ays ten  is 
in:rcas-.d.  Rcl  ;r>..»C'.  y{  *,ivei.  transient  responses  of  the  Buttervortn  filters 
for  first-  throvgh  eighth-order,  are  reproduced  in  Fig.  26. 

W*  ithcr  the  Butte rwort n  :.or  the  binomial  standard  fores  a:v  of  direct  use  in 
Pighv  control  syster  optimization.  However,  they  dc  provide  a  convenient  start¬ 
ing  point  for  devising  other  stai-dard  lores.  The  ITAE  standard  forms  of  Graham 
and  Intlirop  were,  in  fact,  generated  bj  systematically  modifying  Buttervor^h  fores 
of  the  appropriate  order. 
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f’HWTER  IV 


EI7ECT  OK  A  T  LtL\S  TIME  IAG  u?Cn  JIuCCIAL  SiKCfi  KbrwUItH 


This  chapter  presents  t^uertil  exact  for&ulae  expressing  the  effect  of 
system  tire  lags  on  indlcial  error  measures.  These  formulae  cay  be  used  to 
obtain  approximate  relatlors  expressing  performance  measures  for  scoe  fcigh- 
order  systems  In  teres  of  the  performance  of  "equivalent"  lever -crier  systems 
possessing  tic*  lags. 

As  Is  usual,  the  analysis  viil  be  restricted  to  systems  having  zero 
steady-state  error  in  response  to  u  step-input,  to  obviate  infinite  values  of 
indicia!  error  measures.  Three  classes  of  lndlcial  error  measures  will  be 
considered,  designated  aa  its  Kq  16,  l3  of  Chapter  HI,  and  repeated  belov 
for  convenience: 


IU  r 

1  Udl 

(16) 

Jo 

/•CO 

iX'J  - 

J  Wit 

0 

(17) 

r  Ad, 

(18) 

It  vill  be  shovn  that  vhen  z  pure  tin?  d*'»*y,  t,  is  introduced,  the  measures 
can  be  expressed  in  terms  of  the  measures  appropriate  to  the  unddayed  system 


Of  the 

following  formulae: 

IU 

•  V  *  [TOjt-> 

(58) 

ITU 

■•of  1wU  *  H«o 

(») 

rr^u 

■  "of  ^l3Ulx-0  +  *H,.C  ♦  I^U-O 

(6o) 

at 

*  M  T*0 
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Proof: 

IU,  ITU,  IT*TJ,  etc.,  can  be  calculated  by  application  of  the  final  value 
theorem.  Thus, 


TU  • 

3“"=  •  •  S  •  u<‘> 

(61) 

ITU  - 

•  S  '  -  B  ■  u'*> 

(6a) 

T&J  - 

iia  .  1  d2 

—o'  »  d*2"' 

(63) 

For  the  eyaten  vith  In^, 

•J(s)  - 

T  <’  *  •*”  ut^*> 

(6V) 

By  the  fliad  value  theorem,  expanding  e"T8  aa  \  -  is 

1U  . 

/"«it  .  ,^0  .  •  1  •  «(=) 

-  .^To  u<*> 

(65) 

- 

rj.  ♦  lln  u  _(») 

O  a-*o  t»o'  ' 

(66) 

- 

tuo  -  wT-0 

(67) 

vhicn  check*  Eq  >3. 
HU  la  evaluated  aioUarly.  Proa  Eq  62  and  G», 
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Applying  tLc  final  value  theorea 


I  tU{t)dt 


Hr 

s-»o 


11a  U0  j  -t» 

4»— ?r- 


lie 
a  -*o 


(70) 

(71) 


which  checks  Eq  59- 


To  evaluate  If%,  Eq  6}  am  6y  are  ccobined  to  give 

£t%(t)  -  4  w  -  S I  «*”  H  W1  *  «  ”  ¥ 

ds  I  a 

OS 


(«)  - 


-«‘,s 


«w 


.vf  .-«<«> 


(75) 


After  sscc  redaction,  this  yields 

■  4IWU  ’  KJ]t^  *  4  l’o  f  W 

which  checks  Sq  6o. 

EFFECT  CE  TEC  IAS  OS  ITAE  AKD  ITS2  FOR  SKXHD-CKSS  ZSSO-TOSITICS-SSRCR  SYSTQC 

The  effect*  of  a  tine  lag,  c,  on  sooe  performance  aeasures  for  a  second- 
order  lero-pocition-errcr  systen  are  new  considered.  Apart  froo  their  intrinsic 
value,  the  results  also  proride  approximations  to  performance  oeasures  of  nigger- 
order  systems.  The  transfer  functicn  of  a  second-order  systes  with  lag  applied 
directly  to  the  input  is 
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p  2 

8  ♦  2£©  s  +  CO 


(73) 


It  Ij  ccr.vcr.tcit  to  vorit  in  normalized  (nondime ns icnai)  tine  so  that 


c(s)  e~” 

"ft1  *2  ♦  St,  *  1 


De.*,  realization  of  performance  neasu res  is  easily  accomplished  by  use  of 
Eq  26,  27,  tnx  .‘.6  of  Chapter  III- 

For  the  system  of  Eq  7>  with  a  *tep  input 


E(a)  -  r~~ 


Appendix  A  *hcv3  that  for  this  system 


8  +258+1 


.  r. 


xth -  1  -  2(  (78) 


“°  *  _  -=L~co,'\ 

/cost  — -  A  / -  ;2  -  «  rot.  — *L_-  ♦  2  •in*1  A- - -  -— — 

\  2>/l  -  c2  /\  2^'  -  C2  /  2V1  -  t2 


For  {  >  1 


[T«n]T-0  .  [ie]^  .  4 

*  H2  *  1 
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Substitution  of  chc so  expression*  into  Eq  r/j  yields  the  r^A.ijte  graphed  la 
Fig.  22.  Tf*  effect  of  tine  lag  ca  XTE2  has  been  calculated  similarly.  Fr® 
Chapter  HI. 


?hc  effect  of  various  ti&e  delay*  on  ITS2  1*  anoun  la  Fig.  25,  which  La* 
beta  constructed  by  eosbiotrg  Eq  Uh,  ^5,  **>. 

/JFROXiyjlTICB  TO  PEEFORJflUJCE  K5ASUEES  CF  HIGHSt-CRDS* 

CPTXKiU.  2ElC-P0Gm<».J3«QB  SXb'itXS 

Figure  Vj  illustrates  the  ladlcl&l  response  tlac  histories  of  optlsw*  ttav, 
ITS2,  and  IT^E2  r<*ro-pca  it  ion -error  systems.  These  criteria  ray  all  be  regarded 
as  ressesably  valid  and  selective,  at  Jcwil  fur  the  system  oraers  shown.  In 
each  case,  it  will  be  noted  that  the  optlxcra  ntfl -order  systen  response  (where 
n  >  2)  car.  be  fairly  well  approximated  by  adding  a  tics  lag  to  the  opticua 
scccnd-order  response.  This  tine  lag  is  conveniently  chosen  as  equal  to  the 
difference  between  the  delay  tire  of  the  actual  system  and  the  delay  tire  of 
a  zero  tine  lag  secona-order  system.  (Kote  that  the  delay  tine  1*  defined  a* 
tw  for  the  response  to  achieve  50  peiCcul  of  1U  final  value.) 

A  possible  procedure  for  *  yfroxleatlng  to  the  perfursuice  neasures  of 
high-ordcr  systems  now  be  cores  apparent.  The  actual  systen  say  b-  replaced  by 
HU  equivalent  second-order  systen  with  tine  lag,  and  the  perfomance  measure 
calculated  for  the  latter  systen  by  scans  of  Eq  58,  29,  and  60.  This  pro¬ 
cedure  is  In  fact  quite  staple,  and  of  reasonable  accuracy,  as  will  be  demon¬ 
strated  for  the  ITAE  perforsance  measure.  However,  it  does  deiand  knowledge  of 
the  delay  tine,  and  a  netliod  for  estimating  this  pararw'ter  will  now  be  de¬ 
scried. 

ESTIKWIOS  CF  CSXAX  TDC 

Each  of  the  optimal  responses  shown  in  Pig.  15  can  oe  roughly  approxi¬ 
mated  by  a  delayed  nap  function  terminated  near  t:*c  flr*.t  c(t)  rerv  crossing, 
as  sketched  is  Fig.  2*».  If  the  delay  tines  of  the  actual  and  approximated 
response  arc  cade  equal,  ther 


30, 


i'lgure  23.  ITE*  for  ft  Second-Order  Zero-Po*itian- Error  Systea  with  a  Tire  Delay,  t 


?: 


•  T0  +  h cot  7 


(C2> 


Tq  and  y  z*y  be  forties  related  to  the  actual  response  by  de raiding  that  the 
integrated  error  of  tl>e  actual  and  approximate  responses  shall  be  equal.  This 
implies  that 


(T0  ♦  cot  r)  *  T0 
n  - - 5 - 


(8J) 


where  q(  is  the  coefficient  cf  s  in  the  tysten  tiaaafer  function.  (Note:  q,  •  12 
for  a  ur it -nune rater  system)  (See  Chapter  V)  If  the  reap  a  rproxl.  nation  is  valid 
(as  it  Is  for  the  minimum  IT?^,  IT^E^,  IT^E^,  and  ITAE  syatena  of  Fig.  Eq  6} 
vl  11  be  iosely  satisfied,  and  can  be  combined  with  Eq  82  to  give 


*1 


(tt*) 


Figures  2*>  and  26  show  that  Eq  8^  is  closely  satisfied  for  the  nininun  ITAE, 

IT^^j  end  IT^e£  rywtema.  It  also  predicts  the  delay  time  for  sc esc  non  opti¬ 
cal  eystess  with  good  accuracy.  For  example,  Fig-  27  illustrates  the  delay 
tine/q,  relationship  for  the  Futt»rwortn  .  liters  oi  first-  truro^n  elgntr.- order. 
The  lndiclal  responses  of  these  filters  are  graphed  in  Fig.  26.  It  vl  11  be 
observed  tuat  even  lor  tne  elghtn-ordcr  filter,  the  delay  tine  exceeds  qj  by 
only  6.6  percent.  It  is  therefore  conclude*!  lliat  Eq.  8s*  predicts  the  delay  tine 
vitn  good  accuracy  for  optical*  VZKE,  ITE2,  IT2*?,  IT-*?  -jysteaa,  and  for  systecs 
such  as  Puttervorth  filters  which  have  responses  approaching  these  opt ins. 


INTERPRETATION  OF  q, 

The  coefflcle.it  q,  is  equal  U>  the  IE  for  a  unit-nune rator  system,  since 

►  *i 


IE  -  - 

s  -*•  o 


n»1  n-2 

u=  “  *  Vi* 

S  -♦  O  _Ti 


n-» 


q,a  +  1 


tf5> 


qj  (which  is  the  generalized  inverse  velocity  constant  of  Ref.  72)  can  also  be 
interpreted  in  i^rns  of  the  open-loop  systea  characteristics ,  as  follows: 

_ 1  _  1 
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s  ♦  qn.js 
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Figure  25. 


Variation  of  Delay  Tine  with  Integrated  Error 
for  Optisua  ITAE  Zero-Position-Error  System* 


Nondimen  sionol  Time 
Detoy  Time  (Time  To  50%  final  Value) 


Figure  26,  Variation  of  Delay  Tixse  with  Integrated  Error 
for  Optima  IE2,  ITS2,  IT^E2,  and  1T-&2  Syalesa 
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19 


0  5  10 

NONDIMENSIONAL  TIME 

Figure  28.  Indicial  Responses  cf  Butte  rvorth  Fixte.-s  of  First  Through 
Eighth  Orders  (reproduced  fron  Ref.  J?) 
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*r.en, 


sgjjy  *  ♦  0(^2*  ♦  0  +  0, 


(67) 


etc.  are  the  operi-luuy  mots 


I>  cause  for  a  ro realized  scro-pctition-error  ays  ten,  »/G(s)  la  always  of  the 
Torn  s(A}S  ♦  042*  ♦  I)  (!„.:*  ♦  O, 


4,  *  E  (68) 

v.-.erc  K  is  the  : u realized  gain  (or  inverse  velocity  constant). 

K  further  interpretation  can  be  obtained  by  applying  forculae  expressing  the 
coefficients  of  a  polynomial  in  teres  of  its  roots  (e.g.,  Ref.  71*)* 

Let 


sn  *  qr-,sn'’1  Q|8  ♦  1  -  (a  -  oj)(«  -  ag)  (s  -  a^)  (89) 


vfere  a, .  a-,  are  the  closed* lOvy  poles. 

Nov  ala2a3’**a  •  (*On 

and 

°2a3?l*’ * ,an  *  4  ♦  a»a2°3'**cV.i  ■  (*On’  qi 

Dividing  Eq  9»  by  Sq  90 

J_  _l_  J_ 


(90) 


(90 


(92) 


Hence,  for  Oj  ,2  •  •  •  n  *11  real,  is  s  imply  ecjjal  to  the  sue  or  the  sister 
tiae  constants.  (A  second-order  factor  yields  tvo  'tine  constants'  sussing  to 
the  appropriate  2?/a^) 


Sooe  further  approximations  to  delay  tine  are  discussed  later  in  tnls 
chapter. 


Having  obtained  a  good  aoproxlnaticn  to  dMay  tire,  and  interpreted  its 
significance,  the  result  ray  rov  be  applied  to  the  calculation  of  performance 


oC 


neasures. 


EXAM?  12  OF  APPROXIMATE  CALCUIATXCH  OF  FSFttgftNCE 
JC/^mFS  FOR  JSCH-CSCaR  CTSTSJS 


To  the  method  proposed,  it  ;ma  bet*  .JjMef  calculation 

of  ITAS  for  normalized  unit-maeratcr  oy&Urms  of  third-  through  eighth-order. 

The  opticus  response  of  taw;  of  the**  *y*te«s  has  teeu  »pproxlr>ated  by  a  delayed 
second -order  syctea  vltt  a  dasping  ratio  of  0-7-  Figure  29  "hov*  how  the  KA3 
of  such  a  system  varies  witr*  the  tlao  lag.  The  measured  and  predicted  optima 
ITAE  are  caspared  in  Fig.  JO.  The  agreement  is  seen  to  be  good,  especially  for 
syst ess  of  third*  through  sixth-order.  The  ollght  falloff  in  accuracy  obtained 
with  seventh-  and  eighth -order  aysrtecs  i*  u«  to  the  depart  re  of  these  syotexe 
frtf*  the  approximating  fom.  Figure  J1  illustrate*  this  well;  the  egregious 
behavior  of  the  seventh-  and  eighth-order  tine  histories  is  Batched  by  a  cor¬ 
respondingly  unexpected  distribution  cf  poles  (Fig.  21  cf  Ref.  57).  There  dve* 
net  seer:  to  be  any  exploration  of  this  phenomenon.  In  this  connection,  the  de¬ 
parture  cf  the  delay  tiise  free  the  predicted  value  fer  the  seventh-  and  eighth - 
orier  systems  shown  in  Pig.  25  should  also  be  noted.  For  bign-order  systems,  it 
is  important  that  the  effective  t  should  be  estimated  as  closely  as  possible. 
Thus,  vhorcrer  exact  values  of  the  delay  tice  are  available,  they  should  be  used 
in  preference  to  the  approximate  x^  *  ^ .  xais  approximation  is  satisfactory  for 
systems  of  third-  through  sixth-order,  but  underestimates  HAS  for  the  eighth- 
order  aycten  by  20  percent.  Bcvrver,  io  view  cf  thy  large  overshoots  cf  tne 
"optics"  ITAE  seventh-  and  eirhth-ordrr  syotr*«.  it  v**->v  ■-  -*■•* 

discrepant  will  be  significant  in  practical  optimization  problems. 

FURTHER  APPROJHiTIClS  TO  TEE  LEXAY  TIHS 

In  the  course  of  developing  the  eppnxdxaU'*'  described  above,  a  mater  ol 
alternatives  v-re  investigated.  It  is  fait  that  a  brief  discussion  of  two  of 
thtoe  would  be  of  interest. 


The  Measured  jiiaa  values  for  to*  sixth-,  seventh-,  end  eighth- oruer  systems 
were  obtained  by  applying  Sirpsoo'3  rule  and  the  trapezoidal  integration  rule 
tc  the  optima  responses  shewn  in  Fig.  5’,  *nd  by  then  averaging  the  resalts, 
Measured  ITAE  values  for  cpila^  second-  through  fifth-order  systems  are  given 
in  jbf.  57- 


2‘ 


Time  Log  =  Delay  Time- Delay  Time  Far  Second-Order  System 
Damping  Rotio,C*0.7  For  Second-Order  System 


FlCur*  JO.  Coeperiton  of  Xeerured  mod  Predicted  ITAE 


0  6  W  15 


NMKtlMUlMtl  TIM.  *f,t 

Figure  .  IrdiclaX  Eerpccccc  of  ITAE  StcaUrd  Fora*  of  Second*  Through 
Eighth-Order  (reproduced  froa  Ref.  37) 


(A 


for  various 


Pirnt,  delay  tine  was  plotted  against  the  coefficient  of  eQ-1 
optieAl  systems  as  shown  in  Fig.  52.  The  relation  is  not  linear,  except  for  the 
blnoeaial  filters  »Lcrc  the  coefficients  of  s  and  s11"^  are  Identical. 

As  shewn  m  F Jty,  a  fairly  linear  relation  exists  for  the  optima  '2AE 
systems  between  delay  tiae  and  system  order;  this  linearity  is  retained  for 
binomial  filters  although  the  slope  changes.  Linearity  for  the  binomial  filters 
w\xild  be  anticipated  since  the  coefficient  of  s2'*  (»  coefficient  of  s)  is  equal 
to  the  system  order.  Pur  the  optima  HI?,  IT^E?,  o’>d  IT^E?  systems,  the  cor¬ 
relation  with  tysti.'  order  is  nut  as  good. 

It  is  concluded  that  both  uf  these  approxiaat ions  yield  lower  accuracy  than 
the  relaliuuship  expressed  by  £q  6h. 
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Fijurt  JJ.  Correlation  <.[  Do  lay  Tint  vltli  Syatea  Order  for  tlorcalire.1  Unlt-lurarator  Standard  }. 


CHAPTER  V 


PERFCrattHCK  HSASURES  FOR  GQCERAL  DETKKMIHISTTC  IJPUTS 

Ttc  atcp  injx.t  has  become  standardized  si  the  input  required  to  generate 
the  responses  free  which  most  performance  measures  are  obtained.  Indicia} 
7{iwiruu*i  ere  easily  genera  tea  on  an  analog  computer,  ej*l  provide  a  coxvenieat 
neans  to  calculate  the  .espouse  to  •**▼>»  input**  »hvvw>*H  use  of 

Maacl's  integral,  or  equivalent  time-series  procedures.  It  would  be  extremely 
convenient  if  performance  re as urea  appropriate  to  inputs  other  than  steps  could 
be  calculated  in  a  similarly  direct  fashion,  (in  general*  it  is  necessary  first 
to  calculate  the  responae,  and  then  to  appxy  toe  measure  to  the  response) .  in 
this  chapter  it  is  shown  how  for  scae  measures  this  procedure  can  be  "short- 
circuited"  by  calculating  the  response  measure  fro*  Measures  applied  to  the 
input  and  to  the  system  indie lei  or  Impulsive  response.  Thus,  instead  of  ap- 
plying  a  measure  to  a  laplace  transfer**  of  high  order,  it  my  be  calculated  by 
algebraic  operations  on  measures  applied  to  two  lap lace  transforms  of  lower 
order.  It  Is  also  shown  how  certain  performance  measures  of  closed-loop  systems 
as/  similarly  be  calculated  directly  frto  perforaarcc  -.eas’ues  for  the  open- 
loop  system. 

These  procedures  constitute  the  first  step  towards  the  establishment  of  a 
"calculus  of  performs ce  measures,"  i.c.,  a  method  of  expressing  performance  meas¬ 
ures  for  complicated  system  in  terms  of  (core  readily  computed)  measures  for 
simpler  systems.  This  chapter  concludes  with  a  brief  discussion  of  error  coef¬ 
ficients. 

IE  FCR  GENERAL  IHPUTS 

The  response  to  a  unit  impulse,  or  Dirac  5  function,  will  be  referred  to 
here  as  the  "impulsive  admittance."  It  is  also  known  as  the  "weighting  function" 
and  "memory  function."  A  system  having  an  impulsive  admittance  which  is  zero  at 
t  -»•  cd  will  now  be  considered.  By  the  final  value  theorem 


where  is  the  system  transfer  function. 
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The  arcs  enclosed  by  the  impulsive  admittance  and  the  t-sxi*  is 


.  X/1’  “<*>«  -  2$,  •  ■ ;  ■  f& 

I  roy0  *  a(-'  R(s)  -  1 
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The  area  e.tclosed  by  the  input  and  the  t-axis  is,  similarly. 


»(•>  *  l*u 


(96) 


The  area  enclosed  by  the  error  response  to  a  general  input  r(t),  satisfying  the 
condition  that  t^^ao  r(  t)  *  0,  is 
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(97) 


Tills  area  la  liie  integrated  error  response  produced  by  the  actual  inpet.  De¬ 
noting  this  area  by  IE  yields  the  simple  formula 


“  “  “admittance  x  “input 


(98) 


The  implications  of  this  formula  for  system  optimization  using  the  IE  criterion 
are  worth,  of  cement.  The  optimization  procedure  is  assumed  to  consist  of  tne 
adjustment  of  system  parameters  so  that  the  TR  <*  *h*  limit  being 

fixed.  Minimization  of  the  DS^j^ttance  will  thus  result  in  the  minimization 
of  the  IE  for  any  specific  n nit#  input- 

Note  that  the  conditions  that  have  been  imposed  on  the  admittance  and  on 
the  input  result  in  a  finite  (or  zero)  IE. 

* 

As  noted  in  Table  IV,  constraints  must  be  imposed  to  avoid  selecting  £  •  0  for 
TR. 
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Integrated  Error  (1H)  tor  Step  Inputs 

There  is  no  essentia  difficulty  in  extending  the  analysis  to  deal  with 
step  inputs  to  zero- position-,  velocity-,  or  acceleration-error  syst e=s,  hut  the 
straightforward  geooetric  interpretation  of  Eq  9®  is  loot.  For  example,  con¬ 
sider  a  zero-pot ition-err or  system  having  the  transfer  function 
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5i«7 
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*  Vis  '  *  V  ♦  1 


n  _  n-1 


*  *  Vi* 
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(99) 


(100) 


For  a  stop  input 


®  *  »— o 


*  *«n-1s  * 


t  3,a  ♦  1  * 


(101) 


To  apply  Eq  98  here,  it  is  necessary  to  replace  the  actual  «y»t**w  by  a 
substitute  which  has  an  impulsive  admittance  identical  to  the  indicia  response 
of  the  original  system.  For  s  »eh  a  system,  the  transfer  function  is 


4^4- 

«V»»; 


♦  q,« 


*  Vi 


•  a 


(102) 
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11s 


+  Vi8 


*  Vi8 


q  C  +8 


(10?) 


Tbc  V£4  t  is  unity  for  a  Dirac  5  function.  Pence  Eq  93  reduces  to  the  trivial 
form 


IE 


TV  »  1 

edaittance  of  substitute  sjrstca 


(104) 


Integrated  Tigs  Moaent  of  Error  (US) 

With  the  sane  restriction*  cn  the  input  aM  admittance  as  in  the  previous 
sections. 
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'  1  fl'  If?!  dt  ‘ 

11&  1  /  _d_  Efs)  \ 
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(no) 


Hence,  En  106  can  be  interpreted  as 


1TL  »  ^^'a^ttance  X  ^^input  *  ^^adaittance  x  ^input 


Oil) 


An  exanple  of  the  calculation  of  Ute  ITE  by  the  use  of  this  remain  now  follows. 
The  implications  of  Eq  111  with  regard  to  optimisation  are  discussed  subsc- 
ouectly. 

Example  of  the  Calculation  of  ITE  for  a  Secona -Order  Systea 

The  following  exacple  illustrates  the  procedure  for  calculating  HE  using 
Eq  HI,  and  checks  the  result  against  that  obtained  by  direct  calculation. 

The  transfer  function  of  the  systea  considered  is 


E(g)  .  _1.  1 

s  ♦  a  a  ♦  b 


(112) 


The  input  is  assuoed  to  be  described  by  F(t)  -  e”0*  -  e”^*  Evaluating  each 
quantity  on  the  right  side  of  Eq  111 
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Similarly, 
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Slallariy, 

^ Input  * 
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(120) 

Inserting  these  results  in  the  general  formula,  Eq  1 11 ,  yields 
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Conventionally. 

transfora 

,  this  formula  can  he  derived  by  direct  application  of  the  laplace 
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vhlcn  Ic  identical  with  the  result  obtained  froc  the  formula 

GSJERAI.  REIA7KSJS  EETVEEH  CERTAIN  MEASURES  OF 
OPEU-LCC?  AliD  CLOSED- LOO?  SYSTEMS 

5bc  relations  oi  t»e  previous  section  have  been  illicitly  expressed  in 
w:—  v.f  closed-loop  parameters.  Equivalent  ypwr.-loop  expression  arc  now 
derived.  These  expressions  enable  the  closed  loop  IE  end  ITE  to  be  calculated 
without  factoring  (or  even  writing)  the  closed-loop  transfer  function. 

Integrated  Error 

Denoting  the  open-loop  transfer  function  by  G(s),  the  closed-loop  transfer 
function,  with  unity  feedback,  relating  error  to  input  is 
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(127) 


Fro  5  the  preview  section 

^ypen-loop  utu»ivtA.'.c%.  * 


(128) 


Similarly, 


„  Ite  1  >  1 

^closed-loop  admittance  *  s-eo  *  8  *  i'  ♦  G(a)  *  1  ♦  0(0)  (129) 


li;;,  as  snsr-n  is  ttc  iriviouc  icc'.Icn, 
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open  loop  admittance 
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Integrated  Time-Masent  of  Error 

For  &  unity  feedback  closed-loop  systea, 

E(*)[!«G(»8  .  E(»)  (152) 

Differentiating  vith  respect  to  s, 

*<»'  •  £o(c}  *  0  ♦«{*<,  -  £f(«)  (i5j) 

Taxing  toe  licit  as  s  0,  and  waking  use  of  the  results  previously  obtained 
for  IE  and  IS,  Eq  1 35  can  be  rewritten  as 
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free  which 


HE  - 


^^lnput  ^  x  ^^open-locy  admittance 
*  *  ^open-loop  admittance 


o;» 


as  noted  in  Chapter  If,  •US  UlKiUlJUV  lUl  optical  systems  **  snail,  and 
ITAE  »si  IAE  can  be  approximated  viti.  -ooa  accuracy  by  irt  and  1£.  tne 

Units  of  validity  «f  thin  approximation,  the  formulae  developed  in  the  present 
chapter  for  llE  and  2  vith  general  inputs,  and  for  open-loop  and  closed-loop 

fores  of  172  and  IE,  can  also  be  applied  to  ITAE  and  IAE.  Corresponding  re- 
o  o 

let ions nips  fox'  IE* ,  Ijx.  >  nod  exact  formulae  for  ITAE  and  IAE  have  not  yet 
been  obtained. 

The  tiae-vcightcd  integrals  of  the  impulsive  response,  referred  to  as 
K^adnittance  l«*  this  cnapter,  are  sickly  related  to  the  dycmlc  error  cuef- 
fidents  discussed  briefly  in  Chapter  H-  The  relationship  vill  nov  be  demon¬ 
strated  (following  Ref.  11  and  72),  and  its  implications  ctudied. 

Error  Coefficients 


The  generalised  dynamic  error  coefficients  are  defined  in  Ref.  72  as  suc¬ 
cessive  coefficients  of  a  paver  s^dee  expansion  of 


£3  ♦  EjS  Ega' 


036) 


Hence 


vtere  Kq,  Ry,  ...  are  the  dynamic  error  coefficients. 
i!(s)  -  f-c"(=)  •  VVO  «  V2r<s>  *  ('57) 


For  an  impulsive  input  R(s)  »  1. 


E(s)  •  Eq  ▼  Eji  ♦  Sgfc2  y  .  . 

f0  .  -  i  E(S)  - 

o 

J ‘«(‘)dt  ■  5^c 3  •  ;  (•  5  t*Wl)  *  -  Ei 
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-  s^0  =  4-5e(“)  - 

as 


(tM) 


T .  t!»e  vixvr  '-•veffJci^''*s  •’■e  directly  proportional  Lc  aisc-vcighteu 
integrals  of  too  icpiisivc  response. 

For  a  seecrd-crder  vnit-macrator  syctea  having  the  transfer  function 


s“  *  “J5c>n* 

p  2 

S  ♦  2£<&.S  ♦  *i>n 
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Taxing  the  Unit  as  s — ►  0 


*4  .  .a 


SiriJariv.  it  ca**  be  shovu  that 
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None  of  these  error  coefficients  results  in  a  satisfactory  criterion  for  an 
impulsive  input  to  toe  second-order  systca  cotusidered.  As  noted  in  Chapter  H, 
error  coefficients  fail  zt  criteria  sc cause  of  their  inability  to  distinguish 

r/*CD 

c(t)dt J  te(l}dt,  etc. 
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CHAPTIS  VI 


ACC*nu«'Y,  <^2inv ITT,  AMD  POtfEl/BfifflGI  BSXme 


In  Chapter  I,  It  was  stated  that  the  dynamic  perfomenc**  of  a  control 
system  or  element  is  assessed  by  considering  stability,  response  to  desired 
inputs,  response  to  unwanted  input 3,  accuracy,  insensitivity  to  parareter 
chances,  and  pcccr/eatr^  demands.  This  report  Is  principally  concerned  with 
the  response  to  desired  inputs.  As  previously  noted,  the  topic  of  response  to 
undesired  inputs  requires  consideration  cf  ststistically-dcccrib'-i  inputs,  and 
falls  outside  the  scope  of  the  present  report.  Stability  measures  have  been 
summarized  in  Table  1;  the  resaining  characteristics  of  accuracy,  inre~*ittvlty 
to  parameter  changes,  and  pover/energy  demands  are  greatly  dependent  upon  the 
detailed  mechanization  mol  aerulynanic  properties  of  the  particular  flight  con¬ 
trol  cysUm  being  assessed,  and  do  not  lend  themselves  readily  to  generalized 
a todies.  This  chapter  tn^refore  presents  only  a  brief  discussion  of  each  or 
these  liiree  topics  from  a  heuristic  viewpoint. 

Accuracy 

Z k  accuracy  is  ezscutiaiiy  the  suppression  of  error,  it  can  be  studied 

by  the  direct  or  integrated  error  w-usurec  discussed  in  previews  chapters.  Host 
cf  these  measures  have  beer,  related  to  zero-position-error  systems.  Equivalent 
systems  are  normally  of  this  fom,  although  the  linearized  model  of  the  actual 
•ysten  cay  possess  a  small  position  error.  This  steady-state  error  is  explained 
below. 

A  unity  reedback  syst^n  will  be  considered,  with  open-loop  transfer  func¬ 
tion  G(s) . 


For  a  step  input 


£(s) 

Rfi) 


1 

i  ►  gU) 


E(t)  11*  1  .  1  1 

t— ♦CD  *"  S—+Q  *  s  !  ♦  0(s)  *  V  ♦  c(o) 


(1W) 


99 


T 


Thus,  the  system  is  &  true  zero-position-error  system  only  for  C(O)  -*>  ao-  In 
deriving  equivalent  *y*t«s,  this  inaccuracy  is  generally  neglected,  because  for 
a  "good*  system,  the  open-loop  amplitude  ratio  at  lev  frcquer>'*i*»s  is  usually 


residual  error  is  -  0.6b  db  for  the  exact  system;  U  la  error  v*i*  neglected  In 
forming  the  equivalent  system.  Tn  general,  it  Is  believed  teat  this  neglect  Is 
Justified;  however,  circumstances  nay  arise  where  steady-state  errors  cf  this 
magnitude  could  be  significant  factors  in  performance  assessment.  In  such  cases, 
the  system  would  either  be  modified  by  adding  integration  to  the  equalization 
to  reaovc  trim  errors,  or  a  steady-state  error  would  be  permitted,  and  its  magni¬ 
tude  considered  together  with  other  criteria  in  Judging  the  over-all  performance. 


Insensitivity  to  Parameter  Changes 

It  is  desirable  that  the  dynamic  performance  cf  c  given  sy*t^»  ahull  not 
be  degraded  by  ch^ges  in  parameters  occurring  either  as  a  **esult  of  a  change  in 
flight  conditions  or  because  of  discrepancies  between  predicted  end  net usl  com¬ 
ponent  characteristics.  In  flight  control  systems,  the  latter  problem  my  become 
particularly  acute,  because  of  the  well-known  uncertainties  in  derivative  eotima- 
tion.  Important  derivatives  such  as  and  Ip  are  often  made  up  of  components  of 
approximately  equal  magnitude  and  opposite  sign;  the  relative  magnitude  of  the 
over- All  derivative  is  thus  very  sensitive  to  small  changes  in  any  of  its  aero¬ 
dynamic  components. 

Formulation  of  the  performance  measure  bring  employed  in  analytic  terms  is 
the  first  step  towards  mastering  this  *-t — ■’ ,  (as  ILU  report)  this 

process  is  carried  only  to  the  stage  of  defining  the  performance  measure  In  terms 
of  the  transfer  function  pole  and  zero  locations.  The  further  step  required  is 
to  define  these  root  locations  in  terms  of  the  aerodynamic  derivatives,  end  auto 
pilot  gains  and  time  constants.  Reference  k  describes  the  approximate  factoriza¬ 
tion  of  conventional  aircraft  transfer  functions  in  such  literal  tens.  To  apply 
the  procedures  outlined  in  this  report,  it  is  desirable  to  study  (by  the  techniques 
outlined  in  Ref.  k)  the  sensitivity  of  the  equivalent  system  parameters  to  uncer¬ 
tainties  in  the  basic  aerodynamic  and  autopilot  characteristics,  k  good  under¬ 
standing  of  the  effects  or  *«*U  change*  in  perimeters  may  also  be  obtained 
through  the  time  rector  method  of  Ref.  1C  and  21 .  Further  analytes  of  sensitivity 
arc  given  in  Rex.  pp  and  VI,  which  investigate  the  changes  in  the  roots  cf  the 
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characteristic  equation  resulting  *tcc  snail  changes  in  its  coefficient r,  and 
in  Ref.  75>  vhlth  uses  root-locus  technics  tc  study  the  sensitivity  of  the 
closed-loop  roots  to  cpen-loon  parameter  changes. 

Rcwcr/Knergy  Demands 

As  noted  in  Ref.  67,  "if  the  gain  in  a  physical  ays  ten  Is  made  large  enough, 
a  point  Is  reached  at  which  the  peak  acceleration  of  the  output  response  exhibited 
in  tne  linear  ROdel  exceeds  that  which  nay  be  physically  obtained  free  the  pever 
f.ctuator  cf  the  actual  systea.  At  tnls  point,  Ute  linear  nodel  my  cease  to  be  a 
*/alld  basis  for  design.  Either  a  nonlinear  cathcsatical  model  must  be  employed, 
or  the  design  procedure  set  be  uodified  so  that,  although  baaed  cn  linear  theory, 
she  pccsit'lity  of  saturation  is  recognized. r 

It  is  therefore  necessary  to  hare  ccce  check  as  to  whether  or  not  saturation 
is  occurrir^.  Frequently,  this  my  be  zccccplisbed  by  examining  the  anroiitnde 
of  the  p^Oc  overshoot  of  the  output  response.  *Ihe  appropriate  systems  charac¬ 
teristics  graphs  given  in  Chapter  n,  and  in  Ref.  1  *  and  <&,  will  be  found  useful 
for  this  purpose. 

Povcr/cnergy  a*-#ends  nay  also  be  of  Interest  as  direct  performance  measures. 
In  space  vehicles  particularly,  stringent  limitations  oust  be  enforced  on  those 
1  actors.  Under  these  conditions,  dynamic  perfonaance  optimization  Is  achieved  by 
meaict  of  combined  or  constrained  criteria.  Caobired  criteria  art  typically  of 
the  fora 

[iudlcial  error  measure  plus  a  constant  times  total  energy  measure! 

equals  a  eintae.  I 

Constrained  criteria  arc  aerely  single  criteria  of  any  of  the  classes  discussed 
earlier  In  this  report  with  limitations  upon  eaximtai  power,  torque,  or  total 
energy. 

Generalized  assessment  of  combined  criteria  is  very  difficult  because,  in 
the  absence  of  any  specific  application,  selection  of  the  constant  becomes  com¬ 
pletely  arbitrary.  However,  such  criteria  may  be  very  valuable  for  a  given  ap¬ 
plication.  It  is  hoped  that  a  discussion  of  the  use  of  these  criteria  may  be 
included  in  a  subsequent  report. 
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Procedure*  eaployed  for  calculating  certain  indlcial  error  aeaeuree  nay  also 
be  employed  to  calculate  energy  requirement*  fwr  specified  synteae.  ?e.rcrence  5^ 
considers  an  inert  la -wheel  attitude  control  eysteta,  and  shows  how  the  corplex 
convolution  sethod  cnplcjed  to  calculate  TK^  can  be  with  little  ccillia* 

tion  to  calculate  the  energy  expressed  Ui  stal  lilting  the  response  to  a  step  In¬ 
put.  Because  the  tine  ht story  of  the  power  required  ts  a  daaped  sinusoid,  the 
1AE  calculation  sethods  of  the  present  report  could  also  he  tilled  to  tbit  case 
to  calculate  ^|p| dt  (which  equals  the  total  energy  when  no  provision  la  wide  for 
recovering  the  kinetic  energy  of  the  inert  la  wheel) • 


102 


SUMWQf  A5D  COKCJjJSIOKS 


1.  perforw*nc^  »***4  ?55c:^tH  '*r1ter<,\  'ft"  ’Inmr  "?r.!tir.t  ccef- 

flcieul  ays  teas  with  deterministic  inputs  fcav*..  been  investigated,  with  particular 
reference  to  flight  control  system.  The  application  of  performance  measures  has 
been  far>Ht*te<5  by  substituting  for  the  actual  flight  ccutiol  system  an  "equiva¬ 
lent"  low-order  linearised  systea  having  similar  dynamic  characteristics.  This 
equivalent  systen  was  constructed  by  dividing  the  actual  system  transfer  function 
into  regions  of  ir.wcrcet  defined  by 

|0(Jo)j  »  1,  ovci  vblch  1  1 

|3(jo)|  <  1,  over  .m=n  •  |o(j<o)| 

|0(»|  i  I 

Th;  fors  of  ir  bbis  last  regior  ^efjros  tii*  dcrloant  ccdtx  of  tte 

closed-loop  systes  response,  end  can  usually  be  closely  approximated  by  a  eyites 
of  first-,  second-,  or  third  order- 

?.  A  critical  and  exhaustive  survey  cf  current  performance  measures  has 
be -n  cn-vhiew* .  Analytic  for na  f  r  ISAS,  IAS,  etc.,  are  presented,  and  a  n'rtb-*- 
of  errors  ir  previously  published  measures  have  been  corrected-  A  complete  cor¬ 
relation  har  been  given  of  crossover  frequency,  bandvidtu,  phase  margin,  peak 
frcqvcr.cy,  ratification  ratio,  tine-to-peak.  peek  overshoot,  a«d  *>'v- 

second-oider  unit-numerator  systems.  Normalised  presentations  are  used  so  that 
practical  1  ini  tat ions  on  the  time  scale  of  the  response  (e.g.,  due  to  power/ 
inertia  restrictions)  cay  be  taken  into  account  separately.  It  is  concluded  that 
minimum  ITAE  and  mini  cum  TTE2  test  satisfy  the  combined  requirements  of  validity, 
selectivity,  and  case  of  application.  The  ITAE  criterion  yields  smooth  indicial 
responses  havinj  little  overshoot,  but  Its  analytic  description  is  complicated. 

Cf  the  other  indicial  error  aecsures  examined,  minims  IE^  has  simple  analytic 
forms,  but  it  selects  poor  indicxil  responses;  IT^E^  responses  are  as  good  as 
those  selected  by  ITAE,  tut  IT-E^  (a-wi  also  IT^E2)  analytic  expressions  are  too 


complicated  for  general  use.  TTE^  selects  moderately  smooth  aud  well -dumped  re¬ 
sponses  (less  good  than  ITAE),  hut  it  possesses  tractable  analytic  forms.  There¬ 
fore,  I7i  ^s  ruicereoded  for  analytic  investigations,  whereas  ITAE  is  preferred 
for  optimization  a-iag  analog  cctnputers. 

J.  The  integrated  errcr-rtspoasc  ami  integrated  t. Lmc- wcxguu.-d  error  re* 
spense  of  doeed-locp  rystems  to  a  general  deterministic  input  have  been  related 
to  the  corresponding  meas’irec  of  the  response  to  the  impulsive  Input,  which  in 
turn  have  been  expressed  in  terms  of  open-loop  parameters. 
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APPENDIX  A 

AXAirriU  EVALUATION  OF  LAE,  ITAE,  AND  IT^AE  PERFORMANCE  MEASURES 
FCR  UNIT-NUMERATOR  SECOND- ORDER  8X8136 


in 


CALCUIAnOH  OP  THE  JUDICIAL  ERROR  RESPORSE 
OP  A  UIIT-ITHERATOR  SECOKD-ORDeR  SYSTEM 


Tte  transfer  font,  tics  of  the  system  eotiidenal 


£(.i  .  d 

R  ♦  25^s  *  eg 


|  (->  -  .  ff  *  g?°fe* 

R  '  >  *C<V  *  4 


for  *  *tep  input 


E(.)  . 


•“  ♦  25«i,*  ♦  <4 


_ _ 

<*  ♦  t^)2  ♦  (af^rrp)2 


where 


EVAUMnor  op  iae 

Let 


bw 


VTTl 


tin  I  o^l  -"f5  t  +  t 


•] 


(A-l) 


(A-S) 


(A-3) 


C{t)  *  ’?2  -  t,S  t  ♦  t]  (A-A) 

♦  -  Un*’  1  V>  *  t5  -  tin"1  V>  -  C2  -  cos*’  { 


(A-5) 


“  ^==  ^  j  *’**  -  t2  t  +  »J|  At 


■•V  -  c2 
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Zitwi  .  ^=|£“ 1.--1  *.  '  i— i  *| 

|  /t/»  i»*/a>  /»5*/<d  -2»/rt»  #»l»*/a> 

V1  •'O  *  2*/ta  J  s/m  ~  5"/® 

vfecrs  tfc*.  iategra;ji  (vhi  -b  has  been  emitted  fer  brevity)  in 


•st  sin  andt  - 


e  (•*  »lj  arr  -  cocos  o:) 


With  tbe  appropriate  limits,  the  results  of  the  integration  are 

fV/»  /  ,  . 

A*/®  4‘  •  “  '  > 

*W®  /  % 


So 


8(.)  -  £if,  ♦  &-« /-  *  (i  ♦  i)  - 


(*-8) 


But  Jor  Rf  jxs/e&j  >  0, 

1  ♦  2e'«/“  ♦  2e*2«*/«  *  ... 


-  -  rr!^-7s '  f  7  ;-«*£  •  coUi  ©  (A-5) 


*«/<D  —  .  n 

b(->  .  [“  w  -  V/m] 

■sS[-*-]-sa 


(A- 10) 


vuicL  ia»  Uk  Luylacc  transform  of  ]b(t)i.  Bow,  because  |e(t)|  •  IbftjJe*^*1, 

EA(s)  -  B(n-Cafe) 

where  EA(o)  Is  the  Iap^ce  transform  of  the  absolute  error.  Moling  the  substitution. 


*«  ‘  Jfr^t 


*(*  ♦  ?<%)/<= 


r  ..  *<*  ♦  ?=»> 

COth  —7 


"I  c  +  2'eu 

-  1  +  ■  p - "5  (A-11) 

J  (•  ♦  C<%)2  ♦ 


the  1AE  may  now  be  found  through  the  use  of  the  final  value  theon 


«*  •  t -  JZoJnr) 
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Hot  1r£  that 
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EVALUATE*  OF  ITAE 

ITAE  It  defined  at 


ITAE 


f  t|«(t)l4t  . 

*  i^o[-S£A(‘)] 


(A-11.) 


where  Ea(«)  It  th~  Laplace  transform  cf  the  absolute  error  (Eq  A-llh  ^  Ea(s) 
is  calculated  at  follows:  let 


*  ♦  •  a 


Then  Eq  A-ll  becomes 
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Making  tl*  substitution* 
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(A-20) 


Both  the  IAE  and  17  AF  hare  been  calculated  and  tempered  with  the  values  obtained 
by  Graham  and  lathi ©p  (Ref.  57)  (fm  analog  computer  m^chanl  ration)  la  Fig.  17. 
The  analytically  obtained  results  agree  well  with  those  of  Ref.  57  for  IAS  and 
show  minor  discrepancies  cn  STAS. 

As  5  1,  the  ITAE  expression  (Sq  A-20)  becomes  indeterminate,  but  It  c+*. 

be  shewn  thst  In  the  limit  It  reduces  to  («£2  -  1  )/o§  ■  5/a§  for  £  1 . 

Additionally,  the  measures  may  be  written  In  terms  cf  the  rormllzed  measures 
by  application  of  Eq  26  and  27  of  the  main  text  of  this  report. 


*  3  ITAE(?.i) 


(*-«) 


CAICUIAII08  Cf  IT^AE 
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*  — ►o 


d2EA(t) 


(A-25) 


As  shown  previously. 


The  differentiation  nay  be  simplified  by  grouping  the  terms  as  shown  ov  the 
brackets,  and  operating  on  a  grajp  at  a  tine. 


Then 


d£» 

•jy  -  -A(8®  +  E) 


(A- 25) 


and 


d2E*(«) 

d.e 


(A-26) 


The  calculation  of  caci.  group  Is  outlined  on  the  following  rages* 
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The  lest  tens  of  Eq  A-48  Is  the  lA, 


n 


rtz  serialized  IT2AS  tusl  IT^E  are  graphed  in  Pig.  18,  and  canpared  with  the 
values  obtained  oy  Grahaa  ana  IMhrop  (Ref.  1>0  ■  A  scaling  error  of  2  is 
detected  in  the  latter  curve. 
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AFPEKDIX  3 

A  METHOD  FOR  EVALUATING  IAE  AHD  ITAE  FOR  THIRD  -ORDER  SYSTEMS 


A  METHOD  FOR  EVALUATING  IAS  JSD  ITAE  FOR  T'lIPD-QRDFS  SYST3B 

This  appendix  presents  an  am  lytic  procedure  for  evaluating  IAE  and  ITAE 
for  third-order  systems.  The  zeros  of  the  error  tlrt  history  for  tbeae  syst cmn 
are  not  equally  spaced,  and  the  procedure  in  Append!*  A  carrot  be  used.  Iu  the 
vork  p relented  here,  the  error  ti* «  history  la  described  by  a  Fourier-like 
aeries,  the  coefficients  of  yhlch  are  1 1 oe -dependent .  A  similar  technique  vaa 
esgioyed  in  Ref.  20  to  describe  the  output  of  a  linear  full-vave  rectifier  sub¬ 
jected  to  r.  dssged  sinusoidal  i/rput  (accord -order  error  response).  This  appeal ix 
extends  the  procedure  of  Ref.  20  to  third-order  reiponiea  consisting  of  *  d*«fted 
sinusoid  plus  an  o^poiienUal  tern,  anu  snows  now  the  resulting  expression  for 
|e(t)|  nay  be  Integrated  to  give  IAE.  The  procedure  fo*  obtaining  ITAE  is  also 
outlined. 

To  examine  the  accuracy  of  the  method  (i.e.,  the  maker  of  barssnlcs 
required  to  give  acceptable  representations  of  the  actual  functions),  several 
examoles  of  the  calculation  of  |e(t)|  are  given.  It  Is  shewn  that  only  three  or 
four  harmonics  need  be  talen  In  moot  cases  to  achieve  an  accuracy  of  vithln  2  or 
j  percent.  The  convergence  of  the  IAE  seriea  it*  even  more  rapid. 

DETERMINATION  OF  THE  TRANSFER  FUHCTIUN  0?  A  FULL-WAVE  RECTIFIER 

A  full-vave  linear  rectifier  has  a  trausfe»  characteristic  shown  in  Fig.  3^* 


wnerc  x  is  the  input  signal. 
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T*:e  transforms  f+(®)  and  f.(®)  can  be  evaxu&ted  as  follows: 
f*(c)  -  |  x  e*^*dx  -  -h. 


! 


f.<©)  -  f  {-xjc^dx  -  f” X  e-^dx  -  — 

»wn  »o  (“©. 


and  h(x)  in  Eq  B-2  tecoaes 

h(y) 


:-Jco 

vhlcb  can  be  condecsea  to  give 

^♦JOO  y 


£'<jw  »*c<Jro 

Jco  '-e-Jeo  1  * 


-.c^jcd  .  » 

^  jl„ 


CHARACTERIZAIICH  OP  C!E  1KPVT  SICKAL,  x(t) 


(B-5) 


(B-6) 


The  input  signal  is  the  indie  ill  •»rror  response  of  a  normalised  un'.t- 
nunerator  third-order  ays ten 


for  R(e)  ■  l/s 


E(s) 


where 


_ 1 _ 

s^  ♦  b0s2  ♦  ces  +  1 


»2  ♦  V  ♦  °Q  »2  ♦  ba»  +  c„ 

*5‘  V2-0--  ■  <.*r)[<.- «)****] 


(n?  ♦  p2)7  _  , 


bc 

=o 


2a  ♦  7 


1  »  2aj2 


(B-'O 


(B-8) 
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from  which 


e(t)  ■  ae  ♦  kc-0*  els  (gt  +  t) 


(B-9) 


r[(a  -  7)z  *  e2 


*  .  fa'/* 

♦  *  tan*1  §  -  tan*1  ■ 


e(t)  »ay  be  written 


x(t)  «  e(t)  *  A(t)  ♦  V(t)  eoa  0(t) 


-n 


(B-10) 


where  A(t)  •  ae* 

V(t)  .  kc_at 

0<t)  .  *«.♦♦-  f 

CHAEACTERIZATICS  OF  tHE  OtOTW  SIGMI, 

The  output  nay  be  obtained  by  substituting  x(t)(Eq  B-10)  In  Eq  B-6. 


y(t) 


J-^ico 

r  i 

C-J® 


«(A  t  V  cos  0)  -<a(A  ♦  V  cos  0) 


■dm  (B-ll) 


The  exponential  tense  nay  be  expanded  using  the  following  unifonaly  con¬ 
vergent  *w»rie« (Jacobi' -Aw»r  formula.  Ref.  52,  p.  18  and  Ref.  20,  p.  202). 


•*  -  E  «*W  co.  e 


IB-12) 


vhere  £a  -  1  (a  .  1,  1,  },  •••) 

^-©  *  l 

Z9(x)  is  the  nodineu  Bessel  *W  ties  sf  the  first  hind. 
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Mftklz*;  use  of  the  relationship  !*(-*’/  »  (-1  )al(z). 


a  u  co.  ^  r”J®  u(«*)[erf  ♦ 

•v(t)  *  ’ — 


(*-:?) 


(  -  off 

»  «  €V 


b(t) 


j(t)  -  E 


r  001  »*  utpjy*  ♦  (-ov^g 


2«J 


5- JtD 
CO 

y(tj  »  V(t)  Vft,®)  co.  as 


(*-H>) 


(B-15) 


,  ,  ,  J»«)kDC0  ♦  Mftt'Wk 

C{t,»)  -  ^  J  - 1 - J5 - — 


5-Jcd 

Because  {:  ro«  Eq  3-2  and  Eq  B-10) 


y(t)  -  |«(t)l  .  |*(t)  ♦  v(t)  COO  ( 

y(t)  -  v(t)  |b{t)  ♦  co.  e(t)| 


(S-16) 


In  effect  |b(t)  +  co.  S(t)|  ha.  been  expended  in  .  Ffcuner-liie  .erics  with 
tiae  'trying  coefficient.,  i.e.. 


E  E.C(t,a)  co.  aS(t) 
raO 


|h(t)  ♦  co.  6(t)| 


(B-17) 
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EVALUATION  0?  C(t,») 


C(t,Ods  (B-lS) 

O-JCD 


vh ere  c(t,{)  «• 

To  evaluate  the  cocff le  .ecta  C(t,a),  first  consider  the  integral  of  c(t,£) 
around  the  contour  shown  in  Pis* 


Figure  }6.  Contour  of  Integration  for  C(t,a)  Coefficient* 


twining  K,,  X£,  K^.  as  follows: 
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,^-jp 
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Then  as  fi  — ►  a> —  -♦  C(t,a).  Because 
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(B-19) 


*.fce  Bessel  function  I^Cz)  nor  be  approximated  by  t*  for  email  values  of  s. 

The  discussion  is  initially  restricted  to  U*»  case  a  >  a,  when  toe  singu¬ 
larity  of  C(t,0  at  the  origin  vanishes,  eid  C(t,E)  becoaec  analytic  inside,  and 
on  the  contour  of  Fig.  36.  It  therefore  follows  froa  Cauchy's  theorea  that 


K,  ♦  Kg  -  Kj  ♦  K*  -  0 


(B-21) 


The  integrals  Kg  *  vil1  new  be  consi<J<-*-ed.  An  asy^totic  expansion  for 
IB(t)  f©r  larse  values  cf  | z|  is 
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(B-22) 


1J1 


i 
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Ecncc,  for  large  value e  of  0 
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Thus,  »2  0  ud,  •lallftrly,  Bence, 
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where 

-  ■  COS  s  *  -  b^ 

*  B  -  b  ♦Vb2  *  ' 

Applying  L’Hospltal's  rule  to  the  second  terns  or  2  sad  W  it  can  be  shown  that 
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l 1**  ^  , 


(B-J9) 


-V! 

I. 


(be  s)  w  <  1 


(BAO) 


b  >  1 


Because  7.  (for  a  even)  and  V  (for  a  odd)  arc  single -valued  analytic  func¬ 
tions  of  a  for  a  >  0,  the  theory  of  analytic  continuation  (Ref.  16)  permits  the 
restriction  a  >  2  to  be  reaoved.  The  output  signal  can  near  be  found  by  coabin- 
lng  E»i  B-15#  B-37,  aoi  B-33. 


y(t) 


ae"7*  ♦  he*0*  sin  (f>t  ♦  *)  ■  <(l) 
fejjbs  -fr  c|*{('oc  ♦  r)  sin  (Bt  ♦  t)}j  ^ 


|  ~  ^  2  cos  a(Bt  *  *) 

\+  T  £  w  tLn  B^t  *  ♦) 

\  *  a** 


>  b  <  1 


(B-41) 


I** 


where 


b 


x 

V  « 

s  *  sin~lb 

c  »  cos  s  »  V1-*2 

„  BC  COS  B8  +  b  sin  08 

"  »  A 

rafa^  -  1) 

y  ^  nc  sin  no  -  b  coo  ex 

aW4  -  1) 

HUMSWCAL  CHECK  OP  EXPEESSICIf  FOR  |e(t)| 

Equation  B-kl  expresses  |e(t)j  for  a  third-order  system  lndlclal  response 
as  a  Fourier-like  series  with  tin*; -dependent  coefficients.  As  will  be  shown, 

IAE  and  ITAE  can  bo  derived  from  this  series  by  a  straightforward  procedure. 

The  usefulness  or  Eq  B-*1  depends  upon  the  rapidity  of  the  convergence  of 
the  series,  i.e.,  bow  many  harmonics  must  be  included  to  obtain  acceptable  accu¬ 
racy.  Guo  vay  of  cl^ckir^  this  point  would  be  to  evaluate  IAE,  and  cospare  it 
with  the  published  values  in  Eef.  37,  or  with  values  obtained  by  direct  inte gra- 
tl  Gu.  uowever,  as  ahowu  in  Fig.  11  of  Ref.  37,  IAE  is  very  Insensitive  to 
parameter  changes.  It  would  appear  therefore  that,  with  a  small  number  cf 
example*,  «  much  ’-or-  sensitive  snd  thorough  check  can  be  obtained  by  ccsoputing 
|e|  at  selected  Instants.  Integration  produces  a  smoothing  effect  so  the  maker 
of  harmonics  required  to  represent  IAE  accurately  will  be  less  than  the  murker 
required  for  U(t)|. 

Two  responses  sre  considered:  the  ITAE  standard  form  and  a  lightly  damped 
response.  For  the  standard  fora,  the  paraoeters  in  Eq  E-39  are 

a  -  1.2 
k  -  0.72*5 

7  -  0.707 

a  m  0-52 

0-1.07 

♦  •  -0.279  rad  or  -l6  deg 
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and  therefore  e(t),  Eq  B-9,  beccaee 


c(t)  •  l.ae"0-'0'1  *•  0.7<*5e"''‘:521-  sin  (1  .CT/t  -  0.379)  (P-ki>) 

For  the  lightly  damped  system  the  following  parameter  values  were  selected: 

a  *  0.4264 
k  -  0.798 
7  -  1.142 
a  -  0.0542 
2  -  0.935 
♦  *  */2  -  0.7bl1 

and  therefore  e{t) ,  Eq  3-9,  becomes 

e(t)  -  0.4264e*1,1^2t  4  o.790e’O  O5lf2t  sin  (0.935t  */2  -  0.7611)  (B-43) 

The  standard  fora  was  examined  at  t  •  3*20,  4.55,  6.15,  and  6.4l 

nondinetrional  sec,  these  values  being  chosen  for  convenience  of  calculation. 

The  lightly  camped  response  was  examined  at  t  -  0.4105  and  3*66  normalised  -'ec. 
me  results  are  illustrated  in  Fig.  37  through  42,  which  show  that  the  use  of 
only  the  first  five  harmonics  gives  values  within  7  percent  of  the  true  |e(t)| . 
By  taking  6  or  7  harmonics  accuracy  of  within  2  or  3  percent  can  be  obtained. 

EYAUimC*  CP  IAE 

The  integral  of  the  absolute  error  can  be  fourd  from  Eq  B-4i  for  y  (t) 

rco 

jr(t)4t  -  Bo*  j;  0  (B-W>) 

a-  0 

where  j 

H 0  -  f  ’ 

•'o 

oo  r  co 

2d,  -  /  |e(t)|dt 

■ -o  '  •'t, 
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Number  Of  Harmonics 


Figure  37-  Calculated  Vain*  of  j«(t) |  for  Third-Order  Optical  ITA£  System 
vt  Ranker  of  Harwonics  at  t  ■  k.55  Hornalltea  Sec 
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Number  Of  Harmonics 


Figure  J8.  Calculated  Value  of  |e(t)|  for  Third-Order  Optical  ITAE  Systen 
v»  Number  of  Harncaica  at  t  »  >.20  Normalized  Sec 


Porcent  Of  True  |e(t)| 


CJs)  i 

S®’u+I  wafCs+.Oi^+JSB^J 

rw-4- 

t=366 


Number  Of  Harmottics 


Figure  -2.  Cal'.ulftted  7olu«  of  je(t)|  fcr  UjSj-1?  B**y,u  Sy»tes 
v*1  Mir^r  •'.*  ;>naonlc«  t  »  *-^6  noi*ralJi*£  $~c 


1*2 


1 


The  various  tores  of  H0  can  be  evaluated  in  o  etraightforvard  Banner. 
Uo  -  ♦  x’/a  It  [cir.  (t  •*  sic  (5  ♦  7)j 


»  -  ei,  +  i 

The  procedure  fer  Integrating  the  Ua  forrjs  is  core  complicated,  and  will  nov  be 
described  in  detail. 

Evaluation  of  Uq 


Fro:  Sq  B  ti  and  Eq  B  V», 


U0  .  f  £  “(b,  ♦  c)e‘atdt 


t  -  t-r, 


Bote  that  a  Is  not  replaced  t»y  p  in  b,  a,  or  c.  TU  resulting  expressions  con 
be  likened  to  Iaplncc  transforms.  Other  advantages  vill  accrue  vhen  HAE  la 


■  '  U0  ”  A0 


A0  • 

b|  •  e*^ 


“0 

*s 


-  Vi  -  if 

3}  *  atiT'bj 

dr  *  r« .  &''* 

*o  Jo 


{*-*#) 


intcgr*tinf.  tfc  first  tens  by  parts 


i  .  _is_  f c-(i»8«)td. ,  A’/v^ 

*0  L  -<<>♦*'  J0  »♦«  4  "It-X'  j° 


(B-^) 


f°“(1  -  e"^. 

■Jo 


»(£.  i*  ’) 


at  -  — =£-^ - -  (Her.  26)  (B-5C) 


Eq  B — -y  fcvtMR.it 


«  “  ’/2)  3/2} 

- - V~ 


r(i)r{go) 

where  B(x,  *o>  •  *>z  .  is  the  Bet*  function 

tnd  where  z  -  x  +  Jy  la  a  general  complex  variable 

Evaluation  of  U, 


Free  &1  B-l»l  and  S-W 

5V  f00 

U,  -  ~  J  (be  e  a)  sin  (Bt 


♦  Vie  at 


(B-5>) 


(B-52) 


P 


J»0) 

I  (bjC,  ♦  »,)  «lo  ({St  R^Je'^dt 

D 


(B-53) 


i,(%)  -  ^,<=,  ♦  «, 

'•,(?)  ‘  jC.  [«,«]  *  /  0>1C|  ♦  s^c'^dt 
/Wi(p  -  JR)  -  Wi(p  ■ 


—  ■  COS  ® 


/Wi(p  -  j@)  -  Vt(p  -  JS)\  /w  (p  -  jai  4  v.  (r-  *  jp)\ 

\  5  /  *  slB*V  2  '/ 


(B-5M 


♦  -  »  »  w. 


The  express  I  cr.  for  Wj  (p)  is  sialisr  t z  that  for  D ^  the  difference  bcisg  the 
replaceaent  of  p  by  p  -  g  in  ooc  trrm,  and  hy  p  +  g  in  the  other  ter*. 


a  *  '/*>  b(2S£’  3/2} 

r,(p)  -  - g - *— i % - 


2k 


(B-5W 


The  Beta  function  with  one  coaplcx  argument  can  be  written  as  either  an 
infinite  aeries  or  an  infinite  product  (Ref.  52).  The  actual  calculation  of  IAE 
for  specific  paraaeter  values  la  sisplor  with  the  Infinite  product,  but  the  Infi¬ 
nite  series  allows  a  alaple  derivation  of  2TAR,  as  will  be  deaonstrated. 

The  Beta  function  with  one  cccplcx  arg^aent  can  be  written 


®  (-On  8o<Ro  -  ’Hso  '  *)  (So  -  »/ 
eo>  *  £ - — 

n  »0 


Tsrr 


-(tln) 


(B-S6) 


wito  real  and  positive  and  z  complex 


1*5 


r(-')n  60(80  *  ’Heo  -  2)---(bo  *  »)\ 


)<>7 


Ai  »  i£  Kj  «  tao“‘  t  ]T  (ten*' - — - —  -  tan** - £  .  A 

2  ^2  P  cto\  P  *  €  »  2ci  P  ♦  «  ♦  2*f«0 +  njy 


«,  -  5/* 

«0  '  f/® 

the  enluotloo  or  A,  sod  ^  can  be  performed  by  operation*  or.  a  £ode 

«ii Jigres  because  sod  consist  of  an  alternating  act  of  poles  and  zeros  If  33 

la  regarded  as  tbc  replace  variable  a.  0*»ly  those  tent  with  breakpoints  less 
then  a  few  tinea  3  need  be  included,  because  each  pair  of  poles  and  zeros  ore 
qu-te  close  together. 

Evaluation  of  Ug 

Proa  Eq  B-^1  and  BAV 


—  f  Z(t)  cos  p(»  *  t'e-a*it 


(B-62) 


.  U„  -  -UK. 


t  -  c  .  I, 
p  “  a 

2Cj  cos  2s,  ♦  b,  sin 


*1 


W1 


=0 


cos  2(Bt  *  ♦  ♦  BTt)e'I>Idl  (B-6J) 


•••3  ■/■<■-** 


cos  2(Bt  +  S  *  BTjJe’^dT 


(B-<*) 


tkS 


V2(t)  .  (1  -  b3)5/3 

v>  •  /  M  ■  ■-%' 


\  -  5/2 


/«2(p  -  J2?'  *  «(p  *  J5»fl)\  /vjp  -  m)  -  V  fp  * 

2A0  -  \  5  /  «•  a*  -  \ - -35 - -y»ln  2» 


C-Eblnlbg  Eq  B-6j  ui  1-66,  uA  ualna  Eq  B -57,  0,  bocooee 


*»2  “  (-  0"  .2(»2  -  1)(«2  -  a)  -  n) 

2A-  *  it  - 


1  ♦  am)  coo  a*  -  rp  tin  a*\ 
(p  ♦  2gt)2  ♦  ip?  / 


Enluotloo  of  Ua,  M  >  3,  a  odd. 


Pro*  Pq  3-»l  uni  3-M 


T  l  jT 

’-J  I]  \  nfa*  -  1)  / 


(B-69) 


r 


^ — -  -  £  J °°  we(t)  tin  *(0T  4  t  ♦  pr^e^dt 


-3 

Cxld 


V  (t)  ■  SW,  8ln  OL81  -  ty  C03  OB1 

w»(f)  ■  />*«] 

I W  (p  -  .wP)  -  V(j»  ♦  >p) 


m(ia  -  1  )U  CO 

-  s. 


^0  B-  3  )  W  (p  -  >p)  ♦  W  (p  4  >p) 


B  Odu  | 


•  sin  0  1 


cos  ma1  snd  sin  sflj  can  be  expanded  by  the  r  on  owing  xonmiise  (ner. 
r  -  * 

i-2r-l 


roe  Ffi,  • 


»>  3 
a  odd 


n>  3 
»  odd 


£  (-  1)r2^'2r',B(fc  -  r  -  1)!  ,n-2r 

r?0  ^g"tyr - c> 


With  expressions  for  c«s  g#1  snd  sin  w^p)  bereoes 

a  -  1 

•-w  ■  IT 

r  ■  v 

where  /  n-»  \ 


(8-70) 


«): 


(B-T») 


(5-78) 


1>0 


(p)  -  — 


5/?)  -  to. 


U  i&  vbLwlneo.  Dy  coabiniiuj  £q  3-7^-  B-57  *na  3-7C.  The  resulting  e*p4e»»i0u 
1 4  r*ttcr  complicated.  and  Is  not  given  here. 

Srsiu*.tion  of  UB  *>  i,  *  » ven 

7ro=  Eq  B*h  and  B-V» 


U  -  £  /**(•*  "*  • ’ b  I’la  mVo»  P(pt  ♦ 

“  *  \  \  V  -  D  / 


(3-73) 


UB  *  5^  y*°  wa(t)  cos  =(07  -  «’)c'r"dT  (B-7*»> 


»  -  *  ♦  ar. 


-r»fn~  -  l):i  on 


2\> 


z  1 


,'Uf  -  .tes)  ♦  W !(p  ♦  3=0)1 


(•I 


«  (p  -  3=P>  -  «  (r  ♦  3o0)l 


33 


1 


(B-75) 
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where 


-  Z 


coa  m>,  an ^  aui  xi.  can  ce  expanded  by  t:yj  following  fcr^-ulac. 


a  >  5 


IC  (-  OV'^yi  -  r  -  1);=1B- 
r  -  0  r!(m  -  2rJ! 


(3-76*' 


8  In  Jt8 

B  >  2 

m  even 


,  .«I*’  S  (-1)r2a-2r-'(n-r-1):»1,a 
:  *  <*  ”  cl  Z0  ¥VC=  -h  -  TJs 


Because  n  (p;  ic  the  Laplace  tranrfera  cf  w  (tj 
a  D 

/CO 

(bc1  coa  KS^  b^  aln  as^  'e**TdT 


V 

n  even 


(B-77) 


(3-78) 


(B--J9) 


v&ere 

<b_(p)  -  f^O  ’  hf 


jg-Crtl 

2 


a/2*l  < - -g-  -0*g(B-2rJ)  A 


)  t'Ft  >  (-  1)  (a  -  2r)Vl  - 


r 


=7*)  *  <-  >)=/S+:(n  -  Zr)5(£X^jd,  5/a) 
-W  '  - - * - - 


By  coabinlnff  &}  B -IQ.  B-Vf*  and  B-**?  the  expression  for  ’JH  can  be  found. 
Again  the  expression  ia  cosplicattd,  and  vDd  i-ot  V  given  ntre. 


CAIC7T/TI0N  CF  IAi  F<£t  CPSmT*  reilvD-CTOER  8ISTEK 


Tic  leijuirru  for  the  calculation  or  IAJB  for  the  thi  rd -order  sys¬ 
tem  are  gltfer  just  prior  tc  So  B-fc2.  The  values  of  the  ^irs*  few  ccoponents  of 
IaS  are 


*i  -  i  .9s 
U0  «•  0.285 
U}  -  -0.109 


(b-8o) 


The  value  cf  IAF  seeled  o jx  j t  ri*j.  11  cf  Ref.  5?  is  approximately  2.05,  end 
the  calculated  valve  of  ITAE  is  therefore  within  apjr cxi«fc>l>  5  or  6  percent  or 
the  experlmer  tally  obtained  value,  ’roe  Inclusion  of  ours  cccponents  of  IAE  would 
Improve  the  cccuracy. 

The  infinite  product  representation  (Sq  £-6* )  for  tae  Sets  function  was  used 
for  the  numerical  calculation  of  IAE.  Only  the  first  four  faet-ov}  ur-v*  used,  and 
the  accuracy  was  .judged  to  be  acceptable. 

EVAEUATICW  Cl?  ITAE 

ITAE  is  defined  by 

/•co 

ITAE  -  I  t|e(t)|dt  fS-8" 

•'o 


Because  the  expression  used  for  jej  is  In  two  p^rts 

-T|  /*co 

-  j  tedt  +  /  t|e|dt 

Jo 

It  is  convenient  to  define 

JrTj  oo  foo 

*«« »  E  ^  -  I  ‘i«l« 

o  a » 0  J Tj 


(B-82) 


Cs-3j) 


is: 


i.t*  ,  each  M.  corresponds  to  in*  b'^1  component  of  the  contribution  of  tbs  second 
f««tl  ip  r’q  to  TTAL 


Tip.  cxpr-.u>>l^iM>  fur  I-  ■»««!  Liu  ilp  lcrr-i  "re  rath^i  involved  lengthy. 
Ihco.rtcal  checks  will  :x>t  b*  attempted,  »*vt  Jlj  ar-1  Wq  will  be  evaluated  In  lit¬ 
eral  itr  boeau ico  the  calculation  of  tine  t*o  expressions  Is  the  uost  difflcuV- 
part  if  evaluating  ITAF. 

^roji  Eq  D-*  I  and  E-8 3, 


jf  tot  - 


Ha, 


f: 


(B-64) 


The  first  of  these  Integrals  is  quite  simple,  but  the  second  Is  rather  io- 
H,  can  be  rewritten  as  (Integrating  *he  first  tors) 


|  -^  Jj  •  e~>x,(l  -  tTi)|  ♦  k  J  L  (sin  ot  cos  *  ♦  cot.  pt 
•  V  t  sin  (pt  ♦  ije'^dt 


sin  ♦)c*< 


(B-851 


The  -ccond  trra  or  Eq  0-O>  is  tiie  Laplatc  tiwnafora  oi  ♦.  sin  pi  ai>X  t  CC4  pt 
with  a  replacing  the  iaplaec  varl*bl*  o.  The  tnird  ter»  car-  be  vastly  ivalvated 
by  Raking  a  cha*g»  o*  v«T*i*blw  l  *■  t  ♦  T| .  Upon  integrating  the  second  tcra,  £\ 
becomes 


a,  -  ^  -  e"7*1  ( t  -  n,)j  ♦  k?  »m  (*  *  oj) 

-  ke^1  J'  (t  ♦  T.t}(sin  pi  cos  $  ♦  "'o  pt  sin  s)-'; 


(0-66) 


Finally, 


H. 


c'71*  (1 


*  k  jjr  sin 


(t  ♦  A5)  - 


sin  (5  ♦  ♦  y 


Sili  (4  ♦ 


(3-flr) 


aj 


tan 


2te 

i--? 


\ 


Evaluation  of  Sq 

Free  Eq  B-41  and  15-33, 

*0  '  T  j£  t(br  ♦  eje^dt 
let 

t  •  T  ♦  Ti 

?  -  a 

H0  -  \>  /”<»  •  T, ><*,«!  ♦  c^e-P'dt 
**  o 


(B-68) 


(B-3?) 


enc  tens*  In  this  expression  arc  sinilar  to  those  for  UQ  except  for  the 
factor  (:  +  T.)  in  the  integrand.  Since  the  expression  can  be  likened  to  a 
Laplace  transfers 


-  -&< 


T,S(> 


(E-JO) 


Uq  is  giv.-n  by  Eq  B-*l ,  but  the  derivative  with  respect  to  p  offtra  icae 
difficulty.  Representation  of  thn  Beta  function  by  an  infinite  series  (S3  £-56,) 
cl  r*d  tfies  the  aifferentlatioa  considerably  tc.-ausc  p  occurs  only  iu  the  ters 
l/(z  +  n). 

,  t 

The  derivative  with  aspect  to  p  of  1,/#  +  u  i»  Msrely  ♦  a)  vita  a 

cultiplicatire  constant .  Accordingly, 


sc  *  -  Elfr  •»  1  £  <-  ,)1W'so<*o  -  1)(.0  -  a)  ...  (.„  -  n! 

*  “  “1£7P - ***& - TZZ-TTTT-^ - 


-0  •02g(ni)(|J+  1  «  n) 


1  00  {-  1)rft.(a.  -  1)(*(  -  2)  (M,  -  n; 

*g2  n^O  h(r«><£  ♦  a)2 


(B-9D 


The  expressions  rcr  c  >  1,  vtll  be  olcilhr  to  Eq  B-S?  asd  E-90 

\  -  ’§+TlU« 

The  desirability  of  the  seizes  representation  fer  the  Beta  function  la  .tri¬ 
dent.  Bote  that  tie  terns  of  Hj.  should  converge  more  rapidly  than  those  af  IL, 
because  there  is  generally  an  addltioral  factor  of  the  fom  l/(t  ♦  a)  in  etch  cvn- 
sitics  with  respect  to  n. 
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APPBBJIX  C 

raMiB  cr  THE  miffi  rncar  consoLwsrEH  used  id  vbc^trktb 
SEE  EWIMIEC  SISTH1  CCTCKTT 


>57 


DLTAlIb  O'  THE  ftWtt  VUChT  (XWTR01.  SyCT3<  QSZD 
TO  n£*^?7mATE  xK£  23*mi£ff?  SYSTEM  C01CFP1 


i  cf  Chapter  1  S.i  '.c*.  u»in  body  of  the  report  expresses  tit  vp- 
looj.  t  '.’fticr  function  reiati.**  pitch  altitude  to  elevator  deflection  Tor  tee 
fighter  urolcne  totalled  in  Table  III-1  uf  Ref.  2.  The  altitude  le  20,000  ft, 
the  vc.gh;  Iz  ^C.I’OC  lb,  ail  '.r«e  airspeed  ’b  66 o  rt/ce-'  (Mach  Mo.  *  C.oJ*).  The 
&iii*xtuae  transfer  function  is  a*  quoted  le  Kef.  2, 


The  servomotor  plut  amplifier  transfer  function  is  estimated  as 


(C-2) 


Jut  equalization  it  describee  by  Kj  (— j  ♦  1).  Combining  this  vith  the  product 
3i  C-'  and  C-2  yields  the  open-loop  transfer  function  of  tfcc  complete  system. 


